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I. INTRODUCTION 

The PHENIX experiment at the Relativistic Heavy 
Ion Collider (RHIC) at Brookhaven National Labora- 
tory was designed to measure the properties of matter 
at extremely high temperatures and densities. Under 
such conditions, the possibility exists to produce states 
of matter that have not been observed and studied in 
the laboratory. Perhaps the best known of these is the 
quark-gluon plasma (QGP), matter in which the quarks 
are not confined within individual baryons but exist as 
some form of plasma of individual quarks and gluons. It 
should be emphasized that the exact properties of this 
matter are not known and that the characterization of 
the deconfined state, if such a state is produced, will 
form an essential part of the RHIC program. 

One fundamental element of the study of ultrarelativis- 
tic collisions is the characterization of the interaction in 
terms of variables such as the energy produced transverse 
to the beam direction or the number of charged particles. 
These variables are closely related to the collision geome- 
try and are important in understanding global properties 
of the system during the collision. 

This paper describes a systematic study of dEx/dr] 
and dN c h/dr) at mid-rapidity by the PHENIX experiment 
at center-of-mass energies y/s NN — 19.6 GeV, 130 GeV 
and 200 GeV. The centrality dependence of dEx/drj and 
dNch/dr] is characterized by the number of participants, 
determined with a Glauber model, and is studied as func- 
tion of the incident energy. dEx/drj and dN c h/dr] results 
for all four RHIC experiments are included as part of this 
study. The data taken at 19.6 GeV are particularly inter- 
esting because a close comparison to the lower energies of 
the CERN SPS program can be made. Comparisons are 
also made to previous experiments at the Brookhaven 
AGS and CERN SPS at center-of-mass energies of 4.8 
GeV, 8.7 GeV, and 17.2 GeV. Finally, an extensive set of 
collision models describing the E? and N c h distributions 
are compared to the existing data. 



in ■ The PHENIX detector subsystems relevant for the 
physics analysis presented in this paper are listed below. 

Charged particle multiplicity was measured with two 
MWPC layers of the Pad Chambers (PC) called PCI, 
and PC3. These are located in both central arms at 
the radii of 2.5 m and 5.0 m from the beam axis. The 
PCs cover the full central arm acceptance and have an 
efficiency for minimum ionizing particles greater than 
99.5%. The position resolution of PCI was measured 
to be 1.7 mm by 3 mm and twice that for PC3. PCI 
and PC3 can distinguish between two particle tracks if 
they strike the detector with a separation greater than 
4 cm and 8 cm, respectively. 

For the transverse energy measurements, a PbSc sam- 
pling calorimeter (EMCal) 4] from the PHENIX cen- 
tral spectrometers was used. The front face of EM- 
Cal is located 5.1 m from the beam axis. Scintillation 
light produced in the PbSc EMCal towers is read out 
through wavelength shifting fibers which penetrate the 
module. The depth of the PbSc calorimeter is 18 ra- 
diation lengths (Xo) which corresponds to 0.85 nuclear 
interaction lengths. The PbSc calorimeter has an energy 
resolution of 8.1%/ y/E(GeV) ©2.1% for test beam elec- 
trons, with a measured response proportional to the in- 
cident electron energy that is within ±2% over the range 
0.3 < E e < 40.0 GeV j|. 

Two identical Beam-Beam Counters (BBC) |5j consist- 
ing of 64 individual Chcrenkov counters with 2 cm quartz 
glass radiators each cover the full azimuthal angle in the 
pseudorapidity range 3.0 < \-q\ < 3.9. These detectors 
provide a minimum biased (MB) event trigger and tim- 
ing and are also used for event vertex determination. The 
vertex position resolution for central Au + Au events was 
6 mm along the beam axis. 

The Zero Degree Calorimeters (ZDC) are hadronic 
calorimeters located on both sides of the PHENIX detec- 
tor. They cover a rapidity region of 1 77 1 > 6 and measure 
the energy of the spectator neutrons with approximately 
20% energy resolution. 

The BBC and ZDC were used for the centrality deter- 
mination. 



II. PHENIX DETECTOR 

PHENIX is one of four experiments located at 
RHIC lj. The PHENIX detector consists of two cen- 
tral spectrometer arms, designated east and west for 
their location relative to the interaction region, and two 
muon spectrometers, similarly called north and south. 
Each central spectrometer arm covers a rapidity range of 
\r/\ < 0.35 and subtends 90° in azimuth. The muon spec- 
trometers both have full azimuthal coverage with a rapid- 
ity range of -2.2 < 77 < -1.2 (south) and 1.2 < 77 < 2.4 
(north). Additional global detectors are used as input to 
the trigger and for global event characterization such as 
vertex, time of event and centrality determination. A de- 
tailed description of the PHENIX detector can be found 



III. DATA ANALYSIS 

The analysis procedures for the dEr/drj and dN c h/dn 
measured at yfs~^=130 GeV are described in and 
respectively. In this paper the analysis was improved 
including: 

• Inflow and outflow corrections done based on the 
identified particle data, as opposed to HIJING. 

• Corrected trigger efficiency of 92.2+gQ% instead of 
92.0±2±1%. 

• Modified definition of the Et as discussed below. 

The results presented here for y / s^"=130 GeV are con- 
sistent with results previously published. 
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The same data samples with zero magnetic held 
were used for both Et and N c h measurements at each 
beam energy. The analyzed numbers of events are ap- 
proximately 40 x 10 3 , 160 x 10 3 and 270 x 10 3 for 
^/^7=19.6 GeV, 130 GeV and 200 GeV respectively. 

The main steps of the analysis procedure are discussed 
below in connection to the systematic errors associated 
with them. Some additional details can be found in 0, 

Emm- 



A. Et analysis 

The transverse energy (Et) is defined as: 

i 

where 9i is the polar angle. The sum is taken over all 
particles emitted into a fixed solid angle in an event. By 
convention, Ei is taken to be E\ ot — ttim for baryons, 
El ot +mN for antibaryons and E\ ot for all other particles, 
where E\ ot is the total energy of the particle and tojv is 
the nucleon mass . 

The Et measurement presented in this paper was per- 
formed using the PHENIX PbSc EMCal. The EMCal 
absolute energy scale was set using the ir° mass peak re- 
constructed from pairs of EMCal clusters. The value was 
checked against a measurement of the minimum ionizing 
peak for charged particles penetrating along the tower 
axis and the energy/momentum (E/p) peak of identified 
electrons and positrons. The uncertainty in the absolute 
energy scale is 3% in the ^/~s~^ =19.6 GeV data and 1.5% 
in the 130 GeV and 200 GeV data. 

The EMCal acts as a thin but effective hadronic 
calorimeter at mid-rapidity at a collider 0. The mean 
hadron momenta in the EMCal acceptance are approx- 
imately 0.4, 0.55 and 0.9 GcV/c for pion, kaons and 
(anti)protons respectively 0]. Most hadrons stop in the 
EMCal, depositing all their kinetic energy (at pt less 
than 0.35 GeV/c for pions, 0.64 GeV/c for kaons and 
0.94 GeV/c for protons). 

The average EMCal response to the different particle 
species was obtained with a GEANT-based 0] Monte 
Carlo (MC) simulation of the PHENIX detector using 
the HIJING event generator. The HIJING particle 
composition and pt spectra were tuned to the identified 
charged particle spectra and yields in Au + Au collisions 
measured by PHENIX Elat VsTT=200 GeV and 
130 GeV. The NA49 resultsgl [H HU were used for 
EMCal response studies for J~s^ =19.6 GeV data. The 



"deposited" Et emc was about 75% of the total Et "strik- 
ing" the EMCal. This value varied in the ±1.5% range 
for different centralities and beam energies. 

The uncertainty in the EMCal response to hadrons 
gave a 3% error to the total Et- This uncertainty was 
estimated using a comparison between the simulated en- 
ergy deposited by hadrons with different momenta and 
from the test beam data Q. An additional error of 1.3% 
at =19.6 and 200 GeV and 1% at ^/s^ =130 GeV 

comes from the systematic uncertainties in the particle 
composition and momentum distribution. 

Et was computed for each event (Eq. ^) using clus- 
ters with energy greater than 30 MeV composed of adja- 
cent towers with deposited energy of more than 10 MeV 2 . 
The energy losses at the EMCal edges and due to energy 
thresholds, 6% each, were estimated with the absolute 
uncertainty 1.5%. 

The first main issue for the Et measurement is the 
correction for losses for particles which originate within 
the aperture but whose decay products miss the EM- 
Cal (~10%). The second is for the in- flow contribution 
(~24%), which is principally of two types: (1) albedo 
from the magnet poles; (2) particles which originate out- 
side the aperture of the calorimeter but whose decay 
products hit the calorimeter. The in-flow component 
was checked by comparing the MC and the measure- 
ments for events with a vertex just at and inside a pole 
face of the axial central-spectrometer magnet, for which 
the calorimeter aperture was partly shadowed. The esti- 
mated contribution of the in-flow uncertainty to the Et 
uncertainty is 3% 0. 

Since Et measurements are based on the sum of all 
cluster energies in the EMCal, random noise even in 
a small portion of the total number of EMCal towers 
(~15,000 in PbSc) may affect the total energy in the EM- 
Cal, particularly in peripheral collisions. This effect was 
estimated by measuring the total energy in the EMCal 
in very peripheral events with the collision vertex inside 
the magnet poles. In this case, the EMCal is fully shad- 
owed and no energy deposit from beam collisions is ex- 
pected. The estimated contribution was consistent with 
zero. The uncertainty from this effect contributes 3.5% 
systematic error to the Et measurement in the most pe- 
ripheral bin of 45%-50% at ^fs^ =19.6 GeV, 10% in the 
most peripheral bin of 65%-70% at y/s~^ =130 GeV and 
6% to the bin of 65%-70% at ^/s^ =200 GeV. The con- 
tribution to the systematic error for central events was 
negligible. 



1 The definition of Ei in our earlier publication 7] is different for 
antibaryons contribution: Ej ot was used instead of E\ ot + mjr. 
The current definition increases the value of Et by about 4%, 
independent of centrality. 



2 In thresholds of 20 MeV and 3 MeV were applied for the 
cluster and for the tower, respectively. Energy losses due to 
thresholds were properly accounted for in both analyses. 
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B. N c h analysis 

In the absence of a magnetic field, the particle tracks 
are straight lines. The number of tracks in the event was 
determined by combining all hits in PC3 with all hits in 
PCI. The resulting straight lines were projected onto a 
plane containing the beam line and perpendicular to the 
symmetry axis of the PCs. All tracks intersecting the 
plane at a radius less than 25 cm from the event vertex 
were accepted. 95±1% of all real tracks in the event point 
back within this radius. 

The complete set of tracks thus formed contains both 
real tracks and tracks from combinatorial background. 
The latter can be determined using a mixed event tech- 
nique in which each sector in PCI is exchanged with 
its neighbor and the resulting combinatorial background 
measured. 

The average combinatorial background from the mixed 
event analysis was subtracted from the data from the real 
events. Several corrections were subsequently applied. 

A correction of 15.3% accounts for nonsensitive me- 
chanical gaps between the PC sectors, inactive electronic 
readout cards, and dead pads in the PCI and PC3 detec- 
tors. The data were also corrected for the PC efficiency 
for an isolated hit, measured to be 99.5% using cosmic 
rays Q. The combined systematic error from these cor- 
rections was estimated to be 2.5% for single east arm and 
2.3% for both east and west arms. 

Track losses from the finite double hit resolution of 
the PCs depend on the event multiplicity. Losses occur 
in both the direct counting of tracks and in the combi- 
natorial background subtraction. These two effects were 
studied in great detail using Monte Carlo techniques. To 
account for the track losses in the real event sample, a 
correction of 15%, 13% and 6% for the 5% most central 
events was applied at ^s NJV =200, 130, and 19.6 GeV, 
respectively. 

Track loss due to the finite double hit resolution re- 
duces combinatorial background in the real events more 
than in the mixed events. The number of tracks in the 
mixed events must be decreased by 3.6% to account for 
this. The uncertainty in the correction related to the fi- 
nite double hit resolution of the PCs is estimated to be 
3.5% of the number of reconstructed tracks in the most 
central events at y/s NN — 200 GeV. This number was 
deduced from the simulation and cross-checked with an 
artificial 50% increase of the double hit resolution of PCI 
and PC3. 

An additional correction is related to the decays of 
charged particles and feed-down from the decays of neu- 
tral particles as discussed in |S| where it was determined 
using the HIJING event generator. In this paper the mea- 
sured composition of the produced particles at different 
centralities is used at t/s~^ =200 and 130 GeV [l3lll6|. 
The correction related to particle decay varies about ±1% 
over the full range of measured centralities. In mid- 
central events it is — 1 ± 2.9% and +1 ± 2.5% at ^s NN 
=200 GeV and 130 GeV respectively. At the lowest RHIC 



energy the correction is based on NA49 [rH fljl [T^ mea- 
surements at close energy =17.2 GeV and is about 
11±5.7% independent of centrality. The difference be- 
tween y/s^ — 19.6 GeV and 130 GeV arises from the 
decrease of the particle momenta and the width of the 
^-distribution at lower energy which affects the number 
of tracks from decays of the particles coming from ad- 
jacent rapidities. The uncertainty is also larger because 
the correction was based on non-PHENIX data. More 
details on the analysis can be found at 0, ITTl IT^ . 

C. Determination of trigger efficiency and N p . 

The distribution of the number of participants (N p ) in 
Au + Au collisions was determined using a Monte Carlo 
simulation based on the Glauber model. The inelastic 
cross section of p+p collisions used in the Glauber model 
was taken to be 31 mb, 41 mb and 42 mb at ,/s„ ~ = 
19.6 GeV, 130 GeV and 200 GeV respectively Hj and 
varied within ±3 mb in order to get the systematic errors. 
The nuclear density profile p(r) was taken as the Woods- 
Saxon parameterization: 

p(r) = 1/(1 + e (— (2) 

where r n is the nucleus radius and d is a diffuseness pa- 
rameter. Based on the measurements of electron scatter- 
ing from Au nuclei [2l|, r n was set to (6.38±0.27) fm and 
d to (0.54±0.01) fm. 

The BBC detectors are located in a region where the 
number of produced particles is proportional to N p at 
^/s^7 =130 GeV and 200 GeV [22J. By comparing mea- 
sured BBC spectra to simulations, the MB trigger ef- 
ficiency was estimated to be 92.2+ 3 ^% at both ^/T^ 
=200 GeV and 130 GeV, with less than 1% uncertainty 
in the difference between these two energies. 

One can also use the BBC (or ZDC vs. BBC) response 
to define centrality for a given event as a percentile of 
the total geometrical cross section. The BBC ampli- 
tude distribution and ZDC vs BBC signals divided into 
centrality classes are shown in Fig. ^ By matching the 
detector response simulation to the data, N p can be as- 
signed to each centrality class. The results for N p vary by 
less than 0.5% depending on the shape of the cut in the 
ZDC/BBC space and whether the BBC alone was used 
as a centrality measure. The larger error in N p comes 
from model uncertainties and can be parameterized as 
ANp/Np = 0.02 + 3.0/N p . 

At ^/s NN =19.6 GeV, the BBC acceptance partially 
covers the Au nuclei fragmentation region where the re- 
lation between the particle production and N p is not well 
known for peripheral events. This makes the MB trigger 
efficiency model dependent. To avoid this problem, an 
approach based on the Glauber model and the Negative 
Binomial Distribution (NBD) was applied to the data 
from the PHENIX central arm. For the centrality asso- 
ciations, the BBC signal can still be used after applying 
the correction described in the text below. 
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The NBD, written as: 

P(n, n, k) = T{n + k)/(T(k)nl) • (/x/fc) n /(l + p/k) n+k (3) 

represents the number of independent trials n that are 
required to get a number of predetermined successes if the 
average number of successes per trial is [i. The parameter 
k is related to the variance of the distribution by the 
equation (<r//i) 2 = 1/k + By associating n with 

the number of particles produced in the event such that 
n = f(Np), the NBD describes the distribution of hits in 
a detector |23|,|24| produced by a given number of N p . In 
the simplest case when n oc N p , (Ni lit ) = fJ>(N p ). Using 
probability weights for N p from the Glauber model, one 
can construct a distribution of the number of hits in a 
detector. The coefficients /z and k can be obtained by 
fitting the constructed distribution to the experimentally 
measured distribution. 

The number of hits in the PCI detector shown in the 
left panel of Fig.|2|was used to determine the trigger effi- 
ciency. Nhu oc dN c h/dr] can be parameterized as scaling 
with the number of participants Np, where a is between 
1.0 and 1.1 as measured by WA98 at the CERN SPS [H . 
The Glauber/NBD fit to the distribution of the number 
of hits in PCI is shown as the solid line. The fitting 
range is constrained above some number of hits, where 
the trigger efficiency is equal to 1. The efficiency as a 
function of the number of hits in the detector can be 
found by taking the ratio of measured and reconstructed 
distributions. This is shown in the right panel of Fig. El 
Intergated over all Nhits the MB trigger efficiency was 
found to be 81.5 ± 3% at ^fs^ =19.6 GeV. The 1% 
uncertainty due to variation of a from 1.0 to 1.1 was in- 
cluded in the systematic error. An uncertainty on the 
difference between yfs~^ =19.6 GeV and 200 GeV was 
of 1.5%. 

A fraction of events missing in the trigger at all energies 
belongs to the peripheral centrality classes outside the 
centrality range discussed in this paper. 

As a cross check, the same procedure was applied to the 



BBC response at 200 GeV. It was found that the MB trig- 
ger efficiency in Au+Au and d+Au collisions agrees with 
the procedure based on a full simulation within one stan- 
dard deviation of the systematic error. In Au + Au the 
N p in the centrality bins determined using Glauber/NBD 
method agree better than 0.5% to the values used in this 
paper. In d + Au for a single nucleon-nucleon collision 
the MB trigger efficiency was found to be 57% consis- 
tent with 52 ± 7% measured for PHENIX p + p trigger 
efficiency at the same energy using different method |26l . 
Finally, the fraction of expected p + Au collisions in the 
d + Au sample agrees with the fraction of the events in 
which the corresponding ZDC detects the spectator neu- 
tron from the deuteron within better than 1.5%. 

As it was stated above, the BBC detector at ^fs^ 
=19.6 GeV covers a part of the Au nuclei fragmentation 
region and its response is not linear with N p . Also, 
the number of hits in BBC has a strong vertex depen- 
dence mainly due to the fact that the BBC samples differ- 
ent parts of the dN c hj dr\ distribution at different vertices; 
see Fig. [3J The asymmetry of north and south BBC am- 
plitudes in the same event was studied to correct for these 
two effects. Around vertex z = the asymmetry between 
number of hits in north BBC N(z) and south BBC S(z) is 
(N(z) - S(z))/(S(z) + N(zj) cx (rfNch/drfyidNch/drj) 
reflects the slope of the ^-distribution at BBC rapidity. 
In order to use the BBC signal for N p determination the 
observed signals were scaled such that the asymmetry be- 
tween north and south is the same as in the most central 
events where the influence of the fragmentation region is 
negligible. The vertex dependence was also corrected for. 
The results of the correction are shown in Fig. |3J 

The corrected BBC response was used for the central- 
ity determination. Based on both data and Monte-Carlo 
simulation, a systematic error of 2% was added to the 
determination of the centrality classes using the BBC 
correction procedure. 
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FIG. 2: Left panel: Glauber/NBD fit (line) to the distribution of the number of hits in the PCI detector at ^/s NN =19.6 GeV 
(circles). Right panel: MB trigger efficiency as a function of the number of hits. The parameterization is to guide the eye. 
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TABLE I: Summary of systematic errors given in %. When 
the range is given, the first number corresponds to the most 
central bin and the second to the most peripheral bin pre- 
sented in Tables IXTOIXyTl 
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FIG. 3: Average number of hits in BBC north vs. event ver- 
tex at different centralities before the correction (full symbols) 
and after correction (open symbols) at ^/s NN = 19.6 GeV. 



D. Systematic error summary 

Table [I] summarizes the systematic errors discussed in 
this section. The "Energy resp." error for the Et mea- 
surements combines the uncertainties in absolute energy 
scale, hadronic response and energy losses on the EMCal 
edges and from energy thresholds. The resulting error 
for each centrality bin is a quadratic sum of the errors 
listed in the Table. 



IV. RESULTS 
A. PHENIX results 

The distribution of the raw transverse energy, Et emc , 
into the fiducial aperture of two EMCal sectors is shown 



in the left three panels of Fig. 0] for three RHIC energies 
y/s^ = 19.6, 130 and 200 GeV. The lower scale rep- 
resents the fully corrected Et normalized to one unit of 
pseudorapidity and full azimuthal acceptance. The lower 
axis in the plot is not labeled beyond 200 GeV to avoid 
confusion between the true shape of the dEx/dr] distri- 
bution and Et as measured using the limited acceptance 
of two EMCal sectors. 

For the measurements at y/s NN = 19.6 GeV and 
200 GeV, five EMCal sectors (with azimuthal coverage 
A<fi = 112°) were used, while only two sectors (Acj) = 
45°) were available during the PHENIX run at y/s NN 
= 130 GeV. Results obtained with different number of 
sectors at the same energy were consistent within 1.5%. 

The right three panels in Fig. 0] show the number of 
tracks reconstructed in the east arm of the PHENIX 
detector after background subtraction and all correc- 
tions. The lower axis corresponds to measured distri- 
butions normalized to one unit of pseudorapidity and 
full azimuthal acceptance. For a similar reason as for 
the Et measurement, the lower axis is not labeled above 
200 GeV in dN ch /drj. 
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FIG. 4: The distribution of the raw Et in two EMCal sectors (left) and the number of tracks in the east arm of the PHENIX 
detector (right) per MB trigger, measured at three energies. The lower axis corresponds to mid-rapidity values of dEr/dri and 
dNch/drj respectively. Distributions of the four 5% most central bins are also shown in each plot. 



For the N c h measurements at ^Js NN = 130 GeV, only 
the east arm was used, while for the other two energies 
the measurements were made using both PHENIX cen- 
tral arms. The results obtained with two arms at ^/s NN 
=200 GeV and 19.6 GeV are consistent with each other 
within 1.5%. 

The distributions shown in Fig. 0] have a characteristic 
shape with a sharp peak that corresponds to the most 
peripheral events. Missing events caused by the finite 
MB trigger efficiency in peripheral events would make 
this peak even sharper than measured. The plateau in all 
distributions corresponds to mid-central events and the 
fall-off to the most central Au + Au events. The shape of 
the curves in Fig. 0] in the fall-off region is a product of 
the intrinsic fluctuations of the measured quantities and 
the limited acceptance of the detector. 

The distributions for the four most central bins 0%-5% 
to 15%-20% are also shown in each panel. The centroids 
of these distributions were used to calculate the centrality 
dependence of dEx/dr] and dN c h/dr] 3 . The statistical 
uncertainty of all mean values (less than or about 1%) 
determined by the width of the distributions are small 
because of the large size of the event samples. 



The magnitude of dEx/drj and dN c h/dr) at mid- 
rapidity divided by the number of participant pairs as 
a function of N p is shown in Fig. [S] and tabulated in 
Tables IxHIHWI The right three panels show the same 
ratio for dN c h/drj at three RHIC energies. 

The horizontal errors correspond to the uncertainty in 
N p , determined within the framework of the Monte Carlo 
Glauber model. The vertical bars show the full system- 
atic errors of the measurements 4 added quadratically to 
the errors of N p . The lines denote the corridor in which 
the points can be inclined or bent. The statistical errors 
are smaller than the size of the markers. The upper panel 
also shows the results of the two lower panels with open 
markers for comparison. 

An important result from Fig. [S] is an evident con- 
sistency in the behavior of the centrality curves of Et 
shown on the left and N c h shown on the right for all 
measured energies. Both values demonstrate an increase 
from peripheral (65%-70% bin) to the most central events 
by 50%-70% at RHIC energies ^^"=130 GeV and 
200 GeV. For the lowest RHIC energy (^5^7=19.6 GeV) 
this increase is at the level of systematic uncertainties of 
the measurement. One can note that results from PHO- 



3 All plotted and quoted numbers correspond to average values in 
each centrality bin or ratios of those averages. 



4 Here and everywhere errors correspond to one standard devia- 
tion. 
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FIG. 5: dEr/dr] (left) and dN c h/drj (right) divided by the number of participant pairs at three RHIC energies. Errors shown 
with vertical bars are full systematic errors. Lines show the part of the systematic error that allows bending or inclination of 
the points. Horizontal errors denote the uncertainty in determination of N p . 



BOS |43j, show that the total charged particle multiplic- 
ity is proportional to N p while the multiplicity at mid- 
rapidity over N p increases with N p , indicating that the 
pscudorapidity distribution gets more narrow for central 
events. 

The ratios of the dE^/dr] and dN^/dr) per participant 
pair measured at different RHIC energies are shown in 
Fig. Eland tabulated in Table IXVll In these ratios some 
common systematic errors cancel. 

The increase in the Et production between 19.6 GeV 
and 200 GeV (with an average factor of 2.3) is larger 
than for N c h (with average factor of 1.9). This is con- 
sistent with an increase in the particle production per 
participant common to both Et and N c h and a ~20% 
increase in (m<r) of produced particles contributing to 
the Et parameter only. See section |A"T1 and (l6llri|. 

The ratio of 200 GeV/19.6 GeV shows some increase 
from peripheral to central events; however the increase 
is marginally at the level of the systematic errors of the 
measurement. 

The ratio of 200 GeV/ 130 GeV is flat above N p ~ 80 
and is equal to 1.140 ± 0.043 for E T and 1.126 ± 0.036 
for N c h in the most central bin. A rather sharp increase 
between N p —22 and 83 in the ratios of both quantities 
is still at the level of systematic uncertainties. 

The ratio of the transverse energy and charged particle 
multiplicity at mid-rapidity as a function of centrality is 
shown in Fig. [7| for the three energies. The upper plot 



also shows the results displayed in the lower panels for 
comparison. 



The ratio Et/N^ 5 , sometimes called the "Global 
Barometric Observable", triggered considerable discus- 
sion [23,123. It is related to the (my) of the produced 
particles and is observed to be almost independent of 
centrality and incident energy of the collisions within the 
systematic errors of the previous measurements. The 
present paper forges a direct link between the highest 
SPS and lowest RHIC energies, making a more quantita- 
tive study of ET/N c h possible. 



The results presented in Fig. [7| and tabulated in Ta- 
bles IXIIUlXVl show that the centrality dependence of 
ET/N c h is weak and lies within the systematic errors 
plotted with lines. There is a clear increase in Et/N^ 
between ^5^7=19. 6 GeV and 200 GeV. The ^/s^ de- 
pendence of the results is discussed below. 



Ex/Nch is used as a shortcut for (cLEt / drf) / (dN c h / drf) at r) ■■ 
in C.M.S.. 
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FIG. 6: Ratios of dEr/drj (left) and dN c h/dr] (right) measured at different RHIC energies. The errors shown with vertical bars 
are the full systematic errors. Lines show the part of the systematic error that allows bending or inclination of the points. The 
horizontal errors denote the uncertainty in the determination of N v . 
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B. Bjorken Energy Density 



The Bjorken energy density |2£ 
ing 



can be calculated us- 



1 dE T 
Aj_t dy 



(4) 



where r is the formation time and A± is the nuclei trans- 
verse overlap area. 

The transverse overlap area of two colliding nuclei was 
estimated using a Monte Carlo Glauber model A± ~ 



<j x <Jy, where a :J 



and Uy are the widths of x and y po- 



sition distributions of the participating nucleons in the 
transverse plane. The normalization to irR 2 , where R 
is the sum of r n and d parameters in a Woods-Saxon 
parameterization (Eq. |2J) , was done for the most central 
collisions at the impact parameter 6 = 0. 

For the transformation from dEx/drj^^ to 
dEr /dy\ y= o, a scale factor of 1.25 ± 0.05 was used, 
see I A II 

The Bjorken energy density for three RHIC energies is 
plotted in the left panel of Fig. [S] and tabulated in Ta- 
bles Q^nH2^3 For the 5% most central collisions, esj • i~ 
was 2.2 ± 0.2, 4.7 ± 0.5 and 5.4 ± 0.6 GeV/(/m 2 • c) for 
^s NN =19.6, 130 and 200 GeV, respectively. These val- 
ues increase by 2%, 4% and 5%, respectively, for the 
maximal 7Vp Qrt =394, as obtained from extrapolation of 
PHENIX data points. There is a factor of 2.6 increase 
between the "SPS"-like energy {^fs~^ =19.6 GeV) and 
the top RHIC energy y/s^ =200 GeV. The compari- 
son of the only published £b?=3.2 GeV/ fm 3 at SPS for 
head-on collisions [3(j and top RHIC energies, assuming 
the same r=l frn/c, reveals an increase in energy density 
by a factor of only 1.8, which may come from an overes- 
timation in the SPS measurement, as shown in the left 
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panel of Fig. 1131 and discussed in section lA~3l 

Another approach is used by STAR in |3l| for the es- 
timate of the transverse overlap area of the two nuclei 
A± ~ Np^ 3 in Eq. This approach accounts only for 
the common area of colliding nucleons, not nuclei. The 
results are different only in the peripheral bins as shown 
in the right panel of Fig. |S1 For a comparison, the same 
panel shows the result obtained by STAR which agrees 
with PHENIX result within systematic errors, though 
displaying a smaller increase of the energy density with 



C. Comparison to other measurements 

Comparison to the results of other experiments is com- 
plicated by several factors. AGS and SPS data were 
taken in the Laboratory (Lab.) system while the RHIC 
data are in the Center of Mass (C.M.S.) system. Since 
r\ and Et are not boost invariant quantities, the data 
should be converted into the same coordinate system. 
Some experiments provide a complete set of identified 
particle spectra from which information about Et and 
N c h can be deduced. For other experiments, additional 
assumptions are necessary for their published values. Ap- 
pendix A describes how such recalculation was done in 
each particular case. 

The PHENIX results for N c h are compared to the data 
available from the other RHIC experiments. This com- 
parison is shown in the left panels of Fig. 

There is good agreement between the results of 
BRAHMS MA Hp, PHENIX, PHOBOS [HEEI and 
STAR |37ll38j using N p based on a Monte-Carlo Glauber 
model. This agreement is very impressive because all four 
experiments use different apparatuses and techniques to 
measure the charged particle production. The systematic 
errors of all results are uncorrelated, except for those re- 
lated to the same Glauber model which are small. That 
makes it possible to calculate the RHIC average and re- 
duce the systematic uncertainty. The averaged results 
from all four RHIC experiments are plotted in the right 
panel of Fig. [5] and tabulated in Table IXVIII See sec- 
tion for the procedure. 

Figure ITUI compares Et results from the PHENIX and 
STAR [33 experiments. The results are consistent for 
all centralities within systematic errors, though STAR 
cLEt 1 1 dr\ per participant pair has a smaller slope vs. N p 
above ~70 participants, and Et /N c h shown in the lower 
panel is consistent for all N p . 

The RHIC run at y/s^—19.6 GeV allows a con- 
nection between RHIC and SPS data to be made. 
The highest SPS energy of 158A GeV corresponds to 
y/s NN — 17.2 GeV in the C.M.S., making a direct compar- 
ison of RHIC and SPS results possible. This comparison 
is shown in Fig. ^] See sections lA 31 through IA 61 for the 
details of the data compilation. 

Several comments should be made about this compari- 
son. For both measured parameters the PHENIX results 



and the SPS results agree. The WA98 results (see sec- 
tion IA 4|) are systematically higher than the results of 
other experiments, especially for dET/drj. However the 
WA98 data has an additional systematic error common 
to all points shown for the last bin. For N c h the relative 
spread of the SPS results is larger than for the RHIC re- 
sults shown in Fig.El though overall the ^/s NJV =17 .2 GeV 
SPS measurements are consistent with the PHENIX re- 
sult at v / s WJV =19.6 GeV. 

Different SPS and AGS experiments made measure- 
ments at lower energies. The combined data of AGS, 
SPS and RHIC provide a complete picture of the cen- 
trality behavior of Et and N c h as a function of the 
nucleon-nucleon energy. The centrality dependence of 
dN c h/dr] at mid-rapidity measured at v /s JVJV =4.8, 8.7 and 
17.2 GeV by different experiments is shown in Fig. ^| 
See Table IXVIII for the summary of these results and 
sections I A 51 through I A 71 for the details of the data com- 
pilation. 

At the highest SPS energy the averaging procedure is 
the same as for RHIC energies and weighted experimental 
errors are scaled with the reduced x 2_ like factor S (de- 
scribed in section |A~2|) reaching the value of 1.5 at some 
points. For the intermediate SPS energy A /s77i7= 8 ■ 7 GeV, 
two experiments NA45 |4(J and NA50 |4l| reported the 
centrality dependence of dN^/drj at mid-rapidity. The 
discrepancy in the measurements is close to three times 
the quadratic sum of their systematic error. However the 
shapes of the two curves are almost the same. NA49 has 
published results (see section lA~3)l which give one point 
in dN c h/dri at N p =352. This point favors the NA45 re- 
sult 6 . The average centrality curve is produced taking 
into account the shape of the centrality curves reported 
by NA45 and NA50 and the single NA49 point. See sec- 
tion lA 81 for the averaging procedure at y / s^"=8.7 GeV. 
The errors are scaled with the factor S, which reaches a 
value of 2.5 at some points. The AGS results are pre- 
sented with a curve produced from the combined results 
of the E802/E917 experiments (see section lA~7j) . The av- 
eraging procedure in this case is a simple rebinning of the 
data. 

The average SPS centrality dependence at 
^5^7=17. 2 GeV shown in the upper panel in Fig. ^] 
and the average curve of the two RHIC experiments at 
y/s NN = 19.6 GeV shown in the lower panel in Fig. El are 
very similar. Less than 5% increase is expected due to 
the difference in the incident energy between the highest 
SPS and the lowest RHIC energies (see section II V D II 
below). 

The average values presented in Figs. and ^| are 
summarized in Table IXVIII 



The NA57 results a t bo th SPS energies are published without 
systematic errors in E3 . They are currently not considered. 



12 



> 

fl> 4 
O 

iff 




1 * * * * ^ 

A* . a ^ ^ • 200 GeV 

? 6 ^ v "130 GeV 

o 19.6 GeV 



100 



200 



300 



> 4 
O 

iff 

W 



130 GeV 



100 



° PHENIX (A x ~ o x o y ) 
o PHENIX (A x ~ Njf ) 
* STAR (A L ~ Nj;' 3 ) 



200 



300 



FIG. 8: es-, ■ r deduced from the PHENIX data at three RHIC energies (left) and using different estimates of the nuclear 
transverse overlap at x/s NN =130 GeV (right). 




1.5 



20 ° GeV d A * ft "If 

ItT 1 ■phorop; 



PHOBOS 

★ STAR 
♦BRAHMS 



130 GeV 



(I 



19.6 GeV 



100 



200 



~ a4 
m 

#3 
32 



1.5 



200 GeV t 

w 



9 « ^ r 



130 GeV 



19.6 GeV 



© RHIC average 
PHENIX 



300 M 

P 



100 



200 



300 |sj 



FIG. 9: Left panel: dN c h/drj per pair of iV p measured by the four RHIC experiments at different energies. The shaded area is 
the PHOBOS systematic error. Right panel: RHIC average values (including PHENIX) compared to the PHENIX results. 



D. Dependence on the incident nucleon energy. 



The data compilation made in the previous section al- 
lows for a detailed study of the charged particle produc- 
tion in heavy ion reactions at different incident energies 
of colliding nuclei. Although the data on transverse en- 
ergy production is not abundant, a similar comparison 
can be made 0, 0] . 



1. Central Collisions 

FigurelTSlshows the energy dependence for the dEr / dr] 
and dN c h/drj production per pair of participants in the 
most central collisions measured by different experi- 
ments. See sections IA 51 through IA 91 for the details of 
the data compilation. 

The results shown in Fig. ^| are consistent with 
logarithmic scaling as described in 

II IH Use 

of the logarithmic function is phenomenological and 
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- =200 GeV. PHENIX systematic errors 
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tematic scaling error. 



is suggested by the trend of the data in the range of 
available measurements. The agreement of the fits with 
the data in both panels is very good, especially in the 
right panel where the averaged values are used for N p 
=350. The single point of NA49 [3(j is excluded from 
the Et fit (see section IA 3J> . The results of the fit 
dX/d-q = (Q.bNp ■ A)ln(^JsNN / \Z s °nn) arc: 
for E T = 2.35±0.2 GeV and A = 0.73±0.03 GeV 

for N ch x/s% N = 1.48 ± 0.02 GeV and A = 0.74 ± 0.01. 



The parameter y/ s^ N =2.35 GeV obtained from the 
Et fit is slightly above, although within 3er from the 
minimum possible value of ^/s NN —2 x a.m.u. =1.86 GeV. 
The measurement closest to it at ^/s JVN =2.05 GeV done 
by the FOPI experiment allows to estimate the amount 
of dET/drj produced to be 5.0 GeV in the most central 
collisions corresponding to N p =359. Section [A 91 gives 
details of the estimate. This does not disagree with the 
extrapolation of the fit but does indicate that the loga- 
rithmic parameterization requires higher order terms to 
describe how the Et production starts at very low y / s^7. 

The right panel of Fig. ^] shows the logarithmic fit 
to the N c h, data. It agrees well with all dN c h/dr] results 
plotted for N p =350. Unlike that for Et, the fit parame- 
ter \/s° NN for N c h is 1.48±0.02 GeV which is lower than 
the minimum allowed ^s NN . This suggests that above 
2 x a.m.u. the N c h production as a function of ^/s Niv 
should undergo threshold-like behavior, unlike the Et 
production which must approach zero smoothly due to 
energy conservation. 

The FOPI measurement at ^/s JVJV =1.94 GeV and 
2.05 GeV agrees with the extrapolation of the fit at en- 
ergy very close to 2 x a.m.u.. It is an interesting result 



that colliding nuclei with kinetic energies of 0.037 GeV 
and 0.095 GeV per nucleon in the C.M.S. follow the same 
particle production trend as seen at AGS, SPS and RHIC 
energies. 

A fit to the charged particle multiplicity shows a fac- 
tor of 2.2 increase in dN c h/dr) per participant in the most 
central events from the highest energy at the AGS (^/s NJV 
=4.8 GeV) to the highest energy at the SPS {y/s NN 
=17.2 GeV) and a factor of 2.0 from the highest SPS 
energy to the highest RHIC energy (yfs^ =200 GeV). 
Assuming the same behavior extends to the LHC highest 
energy v / s^"=5500 GeV one would expect dN c h/drj = 
(6.1 ± 0.13) ■ (0.5Ay and the increase in particle produc- 
tion from the highest RHIC energy to be ~60% for the 
most central events. With the greater energy, the rapid- 
ity width should increase by ~60% i.e. the total charged 
particle multiplicity at LHC would increase by a factor 
of ~2.6 from the top RHIC energy. 

The ratio of ET/N c h for the most central bin as a 
function of ^/s NN is shown in Fig. Note that the 
line shown in the figure is not the fit to the data points. 
Rather, it is calculated from the fits shown in Fig. ^5] 
The calculation agrees well with the data. 

There are two regions in the plot which can be clearly 
separated. The region from the lowest allowed ^/s NN to 
SPS energy is characterized by a steep increase of the 
ET/Nch ratio with ^/s NN . In this region the increase 
in the incident energy causes an increase in the (tut) of 
the produced particles. The second region starts from 
the SPS energies and continues above. In this region, 
the ET/Nch ratio is very weakly dependent on ^/s NN . 
The incident energy is converted into particle production 
at mid-rapidity rather than into increasing the particle 
(m T ). 

The shape of the ET/N c h curve in the first region is 
governed by the difference in the \J s [ j^ N parameter be- 
tween Et and N c h- In the second region it is dominated 
by the ratio of the A parameters in the fits. This ratio is 
close to 1 GeV. Extrapolating to LHC energies one gets 
a Et /N c h value of (0.92±0.06) GeV. 



2. Centrality shape 

Another interesting question is how the shapes of the 
centrality curves of Et and N c h change with ^/s NN . 

One approach previously used in a number of papers 
is to describe the shape of the centrality dependence as 
a sum of "soft" and "hard" contributions such that the 
"soft" component is proportional to N p and the "hard" 
component to the the number of binary collisions N c : 
A x N p + B x N c . A disadvantage of this approach is 
that the contributions called "soft" and "hard" do not 
necessarily correspond to the physical processes associ- 
ated with these notations. Another approach is to as- 
sume that the production of Et or N c h is proportional 
to Np, although the parameter a does not have any phys- 
ical meaning. 
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FIG. 11: dEr/dr/ (left) and dN c h/dri (right) divided by the number of N p pairs measured by PHENIX at y/s NN =19.Q GeV 
(solid markers) and recalculated from the results of the SPS experiments at the highest energy ^/s JVJV =17.2 GeV (open markers). 
The p + p result of NA49 is marked with an open cross. 
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FIG. 12: dN c h/dri divided by the number of N p pairs mea- 
sured by AGS and SPS experiments and the average taken at 
different energies recalculated in the C.M.S. 

The results of B/A and a obtained from the fits to the 
data at different ^/s NN are summarized in Table [nj Al- 
though the numbers tend to increase with beam energy, 
the values presented in Table ITT1 are consistent with each 
other within the systematic errors. 

The availability of higher quality data would make it 
possible to derive a more conclusive statement about the 
shape of the curves plotted in Figs. l9landlT2l With the 



TABLE II: B/A ratio and parameter a from the fit to the 
data. Errors are calculated assuming a change of the slope of 
centrality curves within the limits of the bending errors for 
PHENIX and full errors for the averaged data (Table PTvTH . 
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GeV 
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200 
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1.16±0.06 
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1.14±0.08 


1.17±0.08 


1.14±0.05 


19.6 


1.13±0.07 


1.09±0.06 


l.lOiO.ll 


17.2 






1.11±0.08 


8.7 






1.06±0.13 


4.8 






1.20±0.24 



present set of data usually limited to N p above 50, a 
large part of the centrality curve is missing or smeared 
by systematic errors. To avoid this, one can compare 
Au + Au collisions to p + p (N p =2) at the same energy. 

Figure [TBI shows dN c h/dr]/(0.5N p ) divided by the pa- 
rameterization plotted in the right panel of Fig. The 
top panel shows the most central events with N p = 350. 
All points are consistent with 1 demonstrating an agree- 
ment of the fit to the data. The points are connected 
with a line for visibility. The middle panel shows results 
for mid-central events with N p — 100 connected with a 
solid line. The dotted line is the same line as in the top 
panel for N p = 350. The points for N p — 100 are lower 
than N p = 350 by a factor of 0.8 — 0.9, over the plotted 
range of incident energies. The lower panel shows p + p 
data corresponding to N p — 2 measured by several ex- 
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FIG. 14: Ratio of Et over N c h for the most central events as 
a function of ^Js NN recalculated into C.M.S.. The line is the 
ratio of two fits shown in Fig. 1131 The band corresponds to 
one standard deviation of the combined error. 



periments. Dotted lines are the same as appear in the 
upper two panels for N p = 350 and 100 and the p + p 
parameterizations are taken from [ZiL l45| . In the range 
of RHIC energies these points are lower by a factor of 
0.65 — 0.75 than the most central events. 

These results indicate that the centrality curves nor- 
malized to the most central collisions have a similar shape 
for all RHIC energies within the errors of available mea- 
surements. 



E. Comparison to models 
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FIG. 15: The three panels show dN c h/dr//(.5N p ) divided by 
the logarithmic parameterization from Fig. 1131 The pan- 
els correspond to N p — 350, 100 and 2 (p + p) from top 
to bottom. Au + Au points are connected with lines also 
shown in lower panels for comparison. The Au + Au data 
is tabulated in Table IXVHI p + p data and parameteriza- 
tions dN/dri = 2.5-0.25Zn(s]vjv)+0.023Zn(s]vjv) 2 (solid line) 
and dN/dri = 0.27fa(sjviv) — 0.32 (dashed line) are taken 
from 0,0. 



A variety of models attempting to describe the behav- 
ior of Et and N c h as a function of centrality at different 
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y/s^7 are available. An updated set of model results 
were collected from several theoretical groups to make a 
comparison as comprehensive as possible. 

Figures 1161 through 1181 show the comparison between 
the existing theoretical models 7 and the data for 19.6, 
130 and 200 GeV. Brief descriptions of the models and 
their main characteristics are given below. 

One of the more commonly used Monte Carlo event 
generators is HIJING flBl liq . This model, like sev- 
eral others, uses pQCD for initial minijet production, 
and the Lund string model |47| for jet fragmentation 
and hadronization. HIJING also includes jet quenching 
and nuclear shadowing. This type of model typically has 
two components, a soft part proportional to N p and a 
hard part proportional to N co u, which partly motivated 
the discussion in section llVD 21 There are also so-called 
saturation models which also rely on pQCD and predict 
that at some fixed scale the gluon and quark phase-space 
density saturates, thus limiting the number of produced 
quarks a nd g luons. An example of this type of model 
is EKRT j4£| , which is referred to as a final state satu- 
ration model. In this paper, comparisons are also made 
to another parton saturation type model, KLN 0> an 
initial state saturation model, and also to models related 
to HIJING, namely Minijet g3 and AMPT [5l|- AMPT 
is a multiphase transport model, and extends HIJING 
by including explicit interactions between initial minijet 
partons and also final state hadronic interactions. Mini- 
jet follows the same two-component model as HIJING 
but also incorporates an energy dependent cut-off scale, 
similar to the saturation models. 

The other models are listed briefly below. SSHM and 
SFM did not have a designated short identifier, so they 
were named somewhat arbitrarily here, based on the 
physics the models incorporate. SSHM (Saturation for 
Semi-Hard Minijet) |52l | is also a two-component model: 
pQCD-based for semi-hard partonic interactions, while 
for the soft particle production it uses the wounded nu- 
cleon model. DSM |53|, the Dual String Model, is ba- 
sically the Dual Parton Model [54|, with the inclusion 
of strings. SFM (String Fusion Model) [5^, is a string 
model which includes hard collisions, collectivity in the 
initial state (string fusion), and rescattering of the pro- 
duced secondaries. Finally, there are the hadronic mod- 
els, LUCIFER 56], a cascade model, with input fixed 
from lower energy data, and LEXUS ;,57j, a Linear Ex- 
trapolation of Ultrarelativistic nucleon-nucleon Scatter- 
ing data to nucleus- nucleus collisions. 

The available model results range from predicting (or 
postdicting) dN c h/dri at one energy to predicting both 
dN ch /dn and dE T /dr] at 19.6, 130 and 200 GeV. The 
models have varying success in reproducing the data. 



7 Models are presented as the best fit by the polynomial of the 
lowest degree which is closer than 1% to any theoretical point 
provided by the authors of the models. The polynomial is plotted 
in the range where points are provided. 



In Fig. it can be seen that KLN is among the 
most successful at describing the dN c h/dr] centrality de- 
pendence for all three energies. However, at ^/s NN 
=19.6 GeV the theoretical curve is steeper than the data. 
This results in a reversed centrality dependence relative 
to the data for the 200 GeV to 19.6 GeV ratio. SSHM 
describes the 130 and 200 GeV data well, for centralities 
above N p ~100, which is the approximate limit of appli- 
cability for this and other saturation models. For the less 
central events, the model values are lower than the data. 
At 19.6 GeV, the model values are significantly higher 
than the data. The saturation model EKRT describes 
the central points at both energies but overshoots the 
more peripheral data points and thus does not reproduce 
the general centrality dependence of the data. For the 
non-saturation models included in this figure, Minijet re- 
produces both the overall scale, as well as the centrality 
and energy dependence of the data rather well, while the 
cascade model LUCIFER describes the central points at 
130 GeV well, but undershoots the less central values at 
this energy. 

Most of the models included in Fig. 1171 provided val- 
ues for all three energies: 19.6, 130 and 200 GeV. SFM is 
in reasonable agreement with the 130 and 200 GeV data, 
but gives much larger values than the data at 19.6 GeV. 
AMPT is in overall good agreement with the data for 
the two higher energies, except for the increasing trend 
in dN c h/drj at the most peripheral events, which is not 
seen in the data. At the lower energy the N c h centrality 
behavior is underestimated. LEXUS rather severely over- 
shoots the data for all energies, indicating that nucleus- 
nucleus effects are not accounted for. The HIJING mod- 
els (version 1.37 and a new version with implemented 
baryon junctions, HIJING B-B) only provide points at 
130 and 200 GeV and are in reasonable agreement with 
the data at those energies, but generally give somewhat 
lower values. The curves shown include quenching and 
shadowing implemented in HIJING. DSM describes 19.6 
GeV reasonably well for all centralities, and the more 
central bins for 130 and 200 GeV, but overpredicts the 
values for semi-central and peripheral events. 

Figure ^] shows the results for the models that pro- 
vide data for both dN c h/dn and dEx/dn. For dEx/dr], 
LEXUS and SFM consistently overshoot the data for all 
energies. In the ratio Et/N^ , LEXUS gives values that 
are too low except at the lowest energy 19.6 GeV. That 
might indicate that the hadronization mechanism allows 
too little energy per particle. The SFM gives values that 
are too large, except for the most peripheral bin, which 
suggests that the particles are assigned transverse masses 
that are too large. The HIJING versions and the related 
AMPT model are in reasonable agreement with the data 
for both dE T /dn and E T /N ch s . 



Note that the HIJING versions available at the time the data 
were collected and used for predictions were in worse agreement 
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Also shown in Fig. ^] are the ratios of results at 
200 GeV to 19.6 GeV, and 200 GeV to 130 GeV, for 
dEx/dr]. These results, especially the comparison of the 
200 GeV to 19.6 GeV data, is intended to make a more 
precise check of the y/s NN dependence of the models. 
SFM fails to describe the 19.6 GeV data and thus can not 
describe the energy dependence probed by these ratios, 
unlike LEXUS which however does not agree well with 
the individual data curves for 19.6, 130 and 200 GeV. 
AMPT and the HIJING versions reproduce the values of 
the ratios well, as expected since they are in reasonable 
agreement with the individual curves. AMPT and Hi- 
ding are also successful in describing the Ex/N c h ratio, 
as illustrated in the lower panels of Fig. El 

To summarize, most models reproduce at least some of 
the data fairly well, but most also fail in describing all the 
data. Since the model results typically are given with- 
out systematic errors, it is not entirely straightforward to 
quantify the level of agreement or disagreement with the 
data. Qualitatively, the models that are most success- 
ful in describing both dEx/dr] and dN c h/dr] in terms of 
the overall trends, both regarding centrality dependence 
and energy dependence, are AMPT, and the HI JING ver- 
sions. KLN and Minijet unfortunately do not give infor- 
mation on dEx/dr] but are successful in describing the 
dN c h/di~i results. The dN c h/drj results thus can either be 
described by the initial state saturation scenario (KLN) 
or by the mini-jet models that need an energy-dependent 
mini-jet cut-off scale as described in |46ll50j to reproduce 
the data. 



V. SUMMARY 

This paper presents a systematic study of the energy 
and centrality dependence of the charged particle multi- 
plicity and transverse energy at mid-rapidity at ^/s NJV = 
19.6, 130 and 200 GeV. 

The yields, divided by the number of participant nu- 
cleons, show a consistent centrality dependence (increase 
from peripheral to central) between dEx / dr] and dN c h/dr] 
for all energies. Furthermore, the increase in the ratio 
E T /N ch from 19.6 GeV to 200 GeV is consistent with 
a 20% increase in (%) with increasing ^/s NJV . The ra- 
tio Ex/Nch shows only a weak centrality dependence at 
RHIC energies. 

For the y/s NN dependence, comparisons were made not 
only among RHIC results but also including data from 
lower energy fixed-target experiments at SPS, AGS and 
SIS. A phenomenological fit, scaling logarithmically with 
y?s NN , describes both dEx/dr] and dN c h/dr) for the most 
central collisions well for all energies. 

Using the fit results, one can separate two regions with 
different particle production mechanisms. The region 
below SPS energy is characterized by a steep increase 
in Et/N^ ~(mj') with y/s NN , whereas for the energies 
above SPS Ex/N c h is found to be weakly dependent on 

\/ S NN ■ 



Within the systematic errors of the measurements the 
shape of the centrality curves of dN c h/ di]/ '(0.5N p ) vs. 
N p were found to be the same in the range of RHIC 
energies and to scale with ln(^/sjvjv)- The same must be 
true for Ex because Ex/N c h has a very weak centrality 
dependence. 

Based on the dEx/dr) measurements, the Bjorken en- 
ergy density estimates were performed and €Bj • t was 
determined to be 5.4 ± 0.6 GeV/(/m 2 • c) at ^/sT^ = 
200 GeV for the most central bin. This is in excess of 
what is believed to be sufficient for a phase transition to 
the new state of matter. The energy density increases by 
about a factor of 2.6 from the top SPS energy to the top 
RHIC energy. 

Finally, a comparison between the RHIC dN c h/drj and 
dEx/dr] data and a collection of models was performed. 
A few models, notably HIJING and AMPT, reproduce 
both dEx/dr] and dN c h/drj rather well for several ener- 
gies. 
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APPENDIX A: RECALCULATION OF THE 
NON-PHENIX EXPERIMENTAL DATA 

Comparisons of dEx/drj and dN c h/drj between differ- 
ent experiments can be made only if the results are pre- 
sented in the same coordinate system since these values 
are not boost invariants. In some cases a full set of identi- 
fied particles measured by one experiment can be recalcu- 
lated into Ex and N c h- Each case that involves handling 
non-PHENIX published data is separately explained in 



Figure ^| shows simulated rapidity distributions for 
Et and N c h hi the C.M.S. and Lab. frames. Plots pre- 
sented here are for illustration purposes only. The in- 
variant distributions which do not change their shape 
under transition from Lab. to C.M.S. are drriT/dy and 
dN^/dy, while all others do. 

In the C.M.S. system, the transition from "77" to "y" 
at mid-rapidity requires a scaling factor between 1.2 and 
1.3. An accurate determination of this coefficient from 
the published data of other experiments is not always 
possible, therefore for SPS and AGS energies a coefficient 
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FIG. 19: Simulated Et (left) and N c h (right) distributions in rapidity and pseudo-rapidity units in C.M.S. and Lab. systems. 



of 1.25 was used. Because of the definition of Et used 
in this paper, dEx/dy « dmr/dy, around mid-rapidity, 
where rriT is a quadratic sum of the particle mass and 
transverse momentum: m-r = \J m 2 + piL. 

In the Lab. system dN c h/dy w dN c h/dr) and 
dEr/dy rs dEx/dr] at maximum rapidity. A 1.04 con- 
version factor was assigned to the transition from "77" to 
u y" in the Lab. system. 

An error of 5% was assigned to any converted value. 
This error also absorbs uncertainties on various assump- 
tions used in the calculations. For example, the contri- 
bution of neutral particles to the total Et is assumed to 
be: 

E?P = {Eji + Ej< )/2. 

H/rp H/rp -f- JZ/rp 

E^ + E^, = E' T + E T (Al) 

2. Averaging procedure 

Average values were calculated for N p =25, 50,... 375. 
The centrality bin corresponding to a given N p can be 
different in different experiments. dN c h/drj per par- 
ticipant and the associated error were deduced by a 
weighted average interpolation from the two nearest val- 
ues of each experiment. The closest value was required 
to be within a proximity of 25 participants from the N p 
value. The error bars are multiplied by the S'-factor, 
where S — y/x 2 /n.d.f. if \ 2 /n.d.f. > 1 or S = 1 oth- 
erwise. See the Particle Data Group reference |23 for 
details. 



TABLE III: Results on identified particle yields in the 
most central events published by NA49 at mid-rapidity at 
% /s JVN =17.2 GeV taken from Fig. 6 in reference |l7j. 



Particle 


7T+ 


7T 


K+ 


K~ 


P 


P 


dN/dy 


167 


165 


32 


15 


33 


5 


error 


10 


10 


4 


5 


1.5 


0.5 



The total yields per participant and number of par- 
ticipants are taken from Fig. 10 in reference [l7| . Using 
shapes of the dN/dy distributions shown in Fig. 7 in ref- 
erence 0| for different centrality bins quantities tabu- 
lated in IIIII and IIVI can be converted into dEx / drj and 
dN c h/dr] per participant pair at mid-rapidity. System- 
atic errors on particle yields are given in Table 1 in the 
same reference. The systematic error on this value is 
not mentioned in the paper, therefore they were taken 
from [19j. The results used in this paper are also given 
in Table Hvl 

For the same and lower ^/s NN , the identified particle 
yields and (mr) were reconstructed using formula (1) 
and Fig. 1 in reference [LSI and also Table II and formula 
(1) and (2) in reference [l!j. The data obtained from 
the Tables and the fits are summarized in Table Ivl Us- 
ing dN/dy and (mj) the values of dEx /drj and dN c h/dr] 
were recalculated in the C.M.S. frame. The accuracy of 
the procedure was verified by the consistency of results 
presented in Table E| and Table 1711 

The single Et point in Fig. ^| is taken from |3(j as 
405 GeV and scaled up by 10%, then divided by pairs of 
N p = 390 as explained in the text. This poin t does not 
agree with the value of Et deduced from [13, d, [Pj • 



3. NA49 

The identified particle yields in the most central events 
at mid-rapidity at ^/s NJV =17.2 GeV are shown in Fig. 6 
in reference [17] . 



4. WA98 

The centrality dependence of the Et , N c h and Et /N r „ 
were read from the plots in Figs. 7 and 14 in reference |2E 



21 



TABLE IV: Total yields of identified particles per partici- 
pant and mean momentum at mid-rapidity in different cen- 
trality bins published by NA49 (P) at v /s]^~=17.2 GeV, from 
Figs. 8 and 10 in reference [r^|. Recalculated values (R) are 
plotted in Fig. El El 



p 

error 



{dJT/dy/Np) 



error 

(Pt + 



P 
P 
P 
P 
P 
P 

P 
P 
P 

P 
P 

P 
P 
P 

P error 

P (p p T ) [GeV/c] 
R dE T /dr]/(.5N p ) 
R error [GeV] 
R dN ch /dr!/(.5N p ] 
R error 



) [GeV/c] 
{pf) [GeV/c] 

dN K+ /dy/N p 
error 

<?# + > [GeV/c] 
dN K /dy/N p 
error 

<pD [GeV/c] 
{p p T ) [GeV/c] 
dN p /dy/N p 



[GeV] 



362 
10 
1.65 
0.08 
0.29 
0.27 

0.27 
0.02 

0.57 

0.15 
0.01 

0.57 
0.77 
0.03 
.005 
0.87 
1.47 
0.11 
1.75 
0.12 



305 242 
15 15 
1.64 1.55 
0.08 0.08 
0.31 0.30 
0.31 0.29 

0.27 0.22 
0.02 0.02 

0.54 0.50 

0.15 0.12 
0.01 0.01 

0.55 0.53 
0.74 0.75 
0.03 0.03 
.005 .005 
0.84 0.80 
1.50 1.35 
0.11 0.10 
1.74 1.62 
0.12 0.11 



189 130 

15 10 

1.48 1.40 

0.08 0.08 

0.30 0.30 

0.30 0.30 

0.19 0.16 

0.02 0.02 

0.49 0.53 

0.10 0.09 

0.01 0.01 

0.51 0.55 

0.73 0.70 

0.03 0.03 

.005 .005 

0.75 0.78 

1.29 1.23 

0.10 0.10 

1.54 1.46 

0.11 0.11 



72 2 
5 

1.42 1.42 
0.08 0.08 
0.29 0.28 
0.30 0.28 

0.14 0.10 
0.02 0.02 

0.55 0.45 

0.07 0.06 
0.01 0.01 

0.51 0.42 

0.65 0.54 

0.02 0.02 

.005 .005 

0.70 0.48 

1.15 1.00 

0.11 0.07 

1.44 1.38 

0.12 0.08 



ublished 
19]. The 



TABLE V: Temperatures of the identified particles 
by NA49 at different y/s NN , as extracted from ^ 
yields are results of the fits of the parameterizations given in 
these publications. 



T [GeV] 
error [GeV] 
dN/dy 




A A d 

.127 .122 .127 

.004 .002 .004 

16. 3.5 0.32 

6.1 0.67 0.23 



T [GeV] 
error [GeV] 
dN/dy 
error 



7T 

.179 
0.01 
132. 
7. 



7T s\ 

.179 .230 

.010 .008 

140. 24.6 

7. 1.2 



=12.4GeV 
K~ 
.217 
.007 



11.7 
0.6 



P 

133 
003 
29. 
6.2 



P 
.120 
.001 
0.7 
0.06 



A 
.133 
.003 
17.5 
4.4 



A 
.120 
.001 
0.8 
0.08 



d 
.133 
.003 
0.85 
0.28 



T [GeV] 
error [GeV] 
dN/dy 




P P 

130 .137 

002 .004 

40. 0.28 

5.8 0.08 



A 
.130 
.002 
17.2 
2.9 



A 
.137 
.004 
0.28 
0.08 



d 
.130 
.002 
1.25 
0.37 



TABLE VI: Recalculated NA49 results, as plotted in 
Figs. El and El 



N p 

dE T /d V /(.5N p ) 
dN ch /d V /(.5N p ] 
Et I N c h [GeV] 



[GeV] 



17.2 
363.±10 
1.50±0.11 
1.86±0.08 
0.81±0.06 



12.4 
352.±10 
1.16±0.09 
1.54±0.07 
0.78±0.06 



8.7 
352. ±10 
0.94±0.07 
1.24±0.06 
0.76±0.06 



and converted to the C.M.S. frame. Results are summa- 
rized in Table EII1 



5. NA45 

The NA45/CERES collaboration did not publish 
results for dN c h/dr) as a function of centrality at 
^Js NN = YJ .2 GeV. The data was taken from Fig. 6.5 in 
reference jH^] and a 10% error was assigned based on 
the analysis procedure. The number of participants was 
taken for the corresponding cross section bin reported by 
the NA50 results |41|. At the lower energy, the results 
were originally published in |59| and then N p was sub- 
sequently corrected, see e.g.: [4(4. The results presented 
in Fig. 4 in reference |40| for charged hadrons h~ and 
(h + — h~) were added together to get dh/dr\ and then con- 
verted to dN c h/dr] in the C.M.S. frame. The published 
and recalculated results are summarized in Table IVIIII 



6. NA50 

Results on N p are taken from Tables 1 and 2 in refer- 
ence ^lj and on multiplicity from Figs. 2 and 4 tabulated 
in captions in reference [4l| . The systematic errors are 
mentioned in the text. There is some discrepancy in the 
results of NA50 and NA45 as shown in Fig. El ln this re- 
spect the comparison made in Table 3 in reference |41| is 
unclear. The results were converted to the C.M.S. frame. 
Recalculated values are given in Table ITXl 



7. E802/E917 

The centrality dependence of tt + , K + yields and (ttit) 
were recalculated from Tables V and VI in reference [6(| • 
Number of participants are taken from Table II in the 
same publication. The results are presented in Table Ixl 

K~ / K + ratio was assigned a value of 0.17 for all cen- 
tralities based on Tables II and III in reference [r3l| . 
This is consistent with results reported in Fig. 6 in refer- 
ence and Fig. 11 in reference |H3]. The proton pro- 
duction reported in Table IV in reference |6(j was com- 
pared to measurements reported in Fig. 2 in reference |r34| 
and Fig. 10 in reference 63] for different centrality bins. 
The results are consistent, p/p ratio was assigned a value 
of 0.0003 based on Fig. 11 in reference A 25% en- 
hancement in n~ /n + ratio for low m-r reported in |r35| 
for the most central bin is not clearly seen in Fig. 11 in 
reference |63j for all centralities. Such an enhancement 
would contribute an additional 8-9% to the total particle 
and transverse energy production. This is less than the 
systematic error on the result and the recalculation er- 
ror, and thus this effect is not considered. The resulting 
values shown in Table EJ are recalculated to mid-rapidity 
in the C.M.S. frame. 



22 



TABLE VII: Published (P) WA98 results at ^5^7=17.2 GeV taken from Figs.7 and 14 in reference |25 
results plotted in Figs. Illj 1121 1131 and 1141 Additional systematic errors are shown in the plots. 



and recalculated (R) 



P N p 

P error 

P dE T /drj/(.5N p ) [GcV] 

P error [GcV] 

R dE T /dr,/(.5N p ) [GcV] 

R error [GcV] 



382 
11 
2.09 
0.12 
2.00 
0.12 



357 

9 
2.06 
0.12 
1.97 
0.11 



311 
7 

2.06 
0.12 
1.97 
0.11 



269 

7 
2.05 
0.12 
1.97 
0.11 



234 

6 
2.03 
0.12 
1.95 
0.12 



201 

5 

2.00 
0.12 
1.92 
0.11 



174 

5 
1.95 
0.11 
1.88 
0.11 



148 
4 
1.96 
0.11 
1.88 
0.11 



128 
4 
1.93 
0.11 
1.85 
0.11 



109 
4 
1.87 
0.11 
1.80 
0.11 



91 

4 
1.84 
0.12 
1.76 
0.12 



75 
3 
1.78 
0.14 
1.71 
0.13 



62 
3 
1.73 
0.15 
1.67 
0.15 



49 
2 
1.70 
0.17 
1.63 
0.17 



30 

11 
1.59 
0.18 
1.53 
0.18 



p 


N p 


381 


355 


310 


268 


231 


199 


171 


145 


124 


105 


87 


72 


58 


46 


36 


27 


20 


13 


9 


p 


error 


10 


9 


8 


7 


6 


6 


6 


5 


5 


5 


4 


5 


4 


3 


3 


3 


2 


2 


12 


p 


dN ch /dr]/(.5N p ) 


2.66 


2.64 


2.62 


2.60 


2.58 


2.54 


2.50 


2.49 


2.43 


2.40 


2.34 


2.32 


2.31 


2.29 


2.20 


2.17 


2.13 


2.17 


2.22 


p 


error 


0.09 


0.09 


0.10 


0.10 


0.10 


0.10 


0.11 


0.11 


0.11 


0.12 


0.14 


0.17 


0.19 


0.21 


0.23 


0.25 


0.26 


0.29 


0.31 


R 


dN ch /dri/(.5N p ) 


2.13 


2.11 


2.10 


2.08 


2.07 


2.03 


2.00 


2.00 


1.94 


1.92 


1.88 


1.85 


1.85 


1.83 


1.76 


1.74 


1.71 


1.74 


1.78 


R 


error 


0.07 


0.08 


0.08 


0.08 


0.08 


0.08 


0.09 


0.09 


0.09 


0.10 


0.12 


0.14 


0.16 


0.16 


0.18 


0.20 


0.21 


0.23 


0.25 


P 


N p 








383 


359 


315 


276 


242 


211 


185 


160 


140 


123 


106 


91 


78 


66 


56 


49 


P 


error 








10 


9 


8 


7 


6 


5 


4 


4 


3 


4 


6 


6 


6 


5 


5 


11 


P 


E T /N ch [GeV] 








0.78 


0.78 


0.79 


0.80 


0.80 


0.81 


0.79 


0.79 


0.81 


0.81 


0.81 


0.80 


0.79 


0.78 


0.77 


0.77 


P 


error [GeV] 








0.06 


0.06 


0.06 


0.06 


0.06 


0.06 


0.06 


0.06 


0.06 


0.06 


0.06 


0.06 


0.06 


0.06 


0.06 


0.06 


R 


E T /N ch [GeV] 








0.94 


0.94 


0.95 


0.96 


0.96 


0.97 


0.95 


0.95 


1.00 


0.98 


0.97 


0.96 


0.94 


0.94 


0.93 


0.93 


R 


error [GeV] 








0.07 


0.07 


0.07 


0.07 


0.07 


0.07 


0.07 


0.07 


0.07 


0.07 


0.07 


0.07 


0.07 


0.07 


0.07 


0.07 



TABLE VIII: Published (P) NA45 results at ^7=17.2 GeV 
taken from Fig. 6.5 in reference |58| and at y /s NN —8.7 GeV 



from Fig. 4 in reference |4C 
ted in Figs. [TT] E21 



and recalculated (R) results plot- 









S NN~- 


L7.2 GeV 








p 


bin 


0-2.3 


2.3-5 


5-8 


8-12 


12-18 


18-23 


23-35 


p 


dN ah /dr] 


420 


350 


300 


250 


210 


170 


125 


R 


N p 


360 


331 


300 


264 


220 


179 


132 


R 


error 


10. 


10 


9 


8 


7 


7 


6 


R 


dN ch /dri/(.5N p ) 


1.87 


1.69 


1.60 


1.51 


1.53 


1.52 


1.52 


R 


error 


0.19 


0.18 


0.17 


0.16 


0.16 


0.16 


0.17 








/ S JVJV = 


=8.7GcV 








P 


JV P 


368 


335 


287 


238 


183 


120 




P 


h~ 


129 


113 


94 


78 


58 


38 




P 


error 


15 


14 


12 


11 


9 


9 




P 


h+ - h" 


52 


46 


39 


30 


22 


15 




P 


error 


12 


8 


8 


7 


6 


6 




R 


dN ch /dr,/(.5N p ) 


1.35 


1.30 


1.27 


1.25 


1.21 


1.21 




R 


error 


0.08 


0.08 


0.08 


0.09 


0.10 


0.12 





TABLE IX: Recalculated NA50 results plotted in 
Figs. [TT] and E21 

,7^777=17.2 GcV 



N p 


354 


294 


246 


205 


173 


120 


error 


12 


10 


8 


8 


8 


11 


dN ch /dr}/(.5N p ) 


1.98 


1.98 


1.94 


1.95 


1.95 


2.10 


error 


0.16 


0.16 


0.16 


0.16 


0.16 


0.17 






\/ S JVJV = ^ 


i.7GcV 








N p 


356 


295 


245 


204 


170 


127 


error 


12 


10 


8 


8 


8 


11 


dN ch /dr)/(.5N p ) 


0.95 


0.90 


0.90 


0.89 


0.90 


0.92 


error 


0.10 


0.09 


0.09 


0.09 


0.09 


0.09 



TABLE X: Centrality dependence of the identified particles 
measured by E802/E866/E917 collaborations. Number of 
participant pairs is published (P) in Table II in reference |r30|. 
7r + and K + values are obtained by extrapolation (E) from 
E802 measurement very close to mid-rapidity. Data is taken 
from Tables V and VI in reference [r30| . Proton data is a com- 
pilation of the results taken from Table IV in reference |6fl ] 
and Fig.2 in reference [o4|. Recalculated values (R) plotted 
in Figs, im and IHl 



p 


N p pairs 


181 


168 


152 


134 


113 


89.5 


62.5 


26.9 


E 


dJV* + /dy 


64.5 


56.8 


47.6 


39.6 


33.3 


25.8 


17.8 


6.89 


E 


error 


3.13 


2.55 


2.75 


1.75 


1.68 


1.37 


0.89 


0.28 


E 


<mj. + ) [GeV] 


.398 


.392 


.387 


.385 


.375 


.365 


.362 


.361 


E 


error 


.013 


.013 


.014 


.012 


.041 


.011 


.010 


.011 


E 


dN K+ /dy 


10.6 


9.28 


8.12 


6.17 


4.91 


3.73 


2.22 


0.74 


E 


error 


0.45 


0.37 


0.32 


0.27 


0.21 


0.16 


0.11 


0.04 


E 


(m£ + ) [GcV] 


.809 


.787 


.774 


.785 


.770 


.740 


.743 


.685 


E 


error 


.034 


.032 


.028 


.028 


.030 


.027 


.025 


.021 


E 


dN p /dy 


62.8 


57.0 


49.4 


43.0 


33.7 


25.2 


16.5 


6.2 


E 


error 


1.7 


1.5 


1.4 


1 


1 


1 


1 


1 


E 


<m*> [GeV] 


1.25 


1.26 


1.24 


1.23 


1.21 


1.19 


1.17 


1.14 


P. 


dE T /dr,/(.5N p ) 


.608 


.580 


.527 


.492 


.460 


.426 


.396 


.335 


R 


error [GeV] 


.146 


.138 


.125 


.116 


.111 


.103 


.098 


.131 


P. 


dN ah /dr,/(,5N p ) 


.903 


.865 


.812 


.773 


.751 


.725 


.699 


.621 


P. 


error 


.135 


.130 


.122 


.116 


.112 


.108 


.108 


.207 


P. 


E T /N ch [GeV] 


.673 


.670 


.649 


.636 


.612 


.588 


.567 


.540 


R 


error [GeV] 


.127 


.123 


.120 


.117 


.115 


.111 


.111 


.112 



TABLE XI: Recalculated values from E802/E917 experiments 
plotted in Figs, till HI 



v^717 [ GcV l 


4.84 


4.27 


3.81 


dE T /drt/(.5N p ) [GcV] 


.579±.087 


.498±.075 


.405±.061 


dN ch /dr)/(.5N p ) 


.851±.128 


.787±.118 


.678±.102 


E T /N ch [GcV] 


.680±.068 


.634±.063 


.598±.060 



For the lower y/s NN the information about particle 
yields and (tot) was extracted for 7r + and K + from Ta- 
ble II and I in reference [6(j respectively, for K~ from 
Table I in reference and for p from Fig. 2 in refer- 
ence |64j . The same assumptions as above were made to 
recalculate values plotted in Figs. 1131 and 1141 The num- 
bers are given in Table fxTI 



8. Averaging procedure at ^/s NN =8.7 GeV 

The averaging procedure is slightly different for this 
curve. First the average results of NA45 and NA50 are 
produced. Then at N p =350, this result is combined with 
the NA49 measurement using the weighted error method. 
A scaling coefficient before and after NA49 averaging is 



TABLE XII: Particle yields measured by FOPI experiment at 
mid-rapidity extracted from Fig. 21 in reference |6q|. 



z 


1 


2 


3 


4 


5-6 


dM/d(cos6) 


43. 


12. 


2. 


0.5 


0.25 


error 


4.3 


1.2 


0.2 


.05 


.025 



calculated. The NA45/NA50 combined result is scaled 
by this factor for all values of N p . 

9. FOPI 

The FOPI results for N ch were calculated for 400A 
MeV based on the data plotted in Fig. 21 in reference [r3e| . 
The points were read at the angle corresponding to the 
mid-rapidity angle (6 = 55°) and then converted to 
dN ch /dr] resulting in 39±4 at ^75^7=2.053 GeV. The 
corresponding number of participants for a 42 mb event 
sample is 359 based on Fig. 8 in reference [6!j. Data for 
150 A MeV were compiled based on the comparison be- 
tween Fig. 13 and Fig. 14 in reference and the used 
definition of rapidity y, resulting in dN c h/ drj = 40±5 at 
= 1.937 GeV. 

The estimate of the Et production at 400A MeV is 
made based on a comparison of the total yields of the 
particles with Z = 1 in |68l and yields of protons and 
deuterons published in 70] . That allowed to determine 
number of all pions at mid-rapidity to be 20.6 and num- 
ber of all hadrons with Z=l 15.2. Assuming that the par- 
ticle temperatures ar e eq ual to T=40 MeV (exact num- 
bers are published in |Ealz3), one can estimate that the 
contribution to Et from pions is + 3/2T and from 
baryons is 3/2T, according to the definition of Et used 
in this paper. Resulting number of 5.0 GeV is a lower 
limit estimate because the contribution of heavier par- 
ticles is not considered. A conservative error of 30% is 
assigned to this number. 

10. PHOBOS measurement at ^Ts^ =56 GeV 

The PHOBOS experiment published dN ch /drj 
=408±12(stat)±30(syst) at ^/s^ = 56 GeV measured 
for N p =330±4(stat)±}°(syst) in In the same paper 
dN c h/dri per participant between 130 GeV and 56 GeV 
is measured to increase by 1.31±0.04(stat)±0.05(syst). 
That allows for use of the averaged value at ^s NN 
=130 GeV consistent with the PHOBOS result pub- 
lished in [35| to recalculate dN c h/dri at y/s NN = 56 GeV 
with smaller systematic error. This value is plotted in 

Fig. ma 



APPENDIX B: OUTPUT TABLES. 
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TABLE XIII: Results of the measurements by PHENIX at ^/s NN =200 GeV. Errors have the same dimension as preceding 
value. Results are plotted in Figs. |S] |7| and |H] 



bin [%] 




0-5 


5-10 


10-15 


15-20 


20-25 


25-30 


30-35 


35-40 


40-45 


45-50 


50-55 


55-60 


60-65 


65-70 


N p 




353. 


300. 


254. 


215. 


181. 


151. 


125. 


103. 


83.3 


66.7 


52.5 


40.2 


30.2 


22.0 


syst. 


error 


10. 


9.0 


8.1 


7.3 


6.6 


6.0 


5.5 


5.1 


4.7 


4.3 


4.1 


3.8 


3.6 


3.4 


A ± [fm 2 ] 




140. 


125. 


112. 


100. 


90.8 


82.2 


73.9 


66.8 


60.0 


54.3 


49.3 


45.1 


40.9 


37.5 


syst. 


error 


11. 


10. 


9.1 


8.2 


7.4 


6.8 


6.2 


5.7 


5.2 


4.8 


4.6 


4.5 


4.5 


4.9 


dE T /dr] [GcV] 




606. 


493. 


402. 


328. 


266. 


216. 


173. 


137. 


107. 


81.8 


60.4 


43.9 


31.1 


21.1 


stat. 


error 


0.6 


0.6 


0.5 


0.5 


0.4 


0.4 


3 


3 


0.2 


0.2 


0.2 


0.1 


0.1 


0.1 


bending syst. 


error 


2.4 


5.1 


7.0 


7.6 


8.2 


8.5 


8.2 


8.1 


7.6 


6.9 


6.4 


5.5 


4.6 


3.9 


full syst. 


error 


32. 


27. 


22. 


19. 


16. 


14. 


12. 


11. 


9.5 


8.1 


7.2 


5.9 


4.9 


4.0 


e Sj r [GcV fm~ 2 . 




5.4 


4.9 


4.5 


4.1 


3.7 


3.3 


2.9 


2.6 


2.2 


1.9 


1.5 


1.2 


1.0 


0.7 


full syst. 


error 


0.6 


0.5 


0.5 


0.5 


0.4 


0.4 


0.4 


0.3 


0.3 


0.3 


0.2 


0.2 


0.2 


0.2 


dE T /dr]/{0.bN p ) IGcVj 


3.43 


3.28 


3.16 


3.05 


2.94 


2.86 


2.76 


2.66 


2.57 


2.45 


2.30 


2.18 


2.06 


1.92 


bending syst. 


error 


0.01 


0.03 


0.05 


0.07 


0.09 


0.12 


0.14 


0.17 


0.19 


0.22 


0.27 


0.31 


0.36 


0.43 


full syst. 


error 


0.21 


0.20 


0.20 


0.20 


0.21 


0.22 


0.23 


0.25 


0.27 


0.29 


0.33 


0.36 


0.41 


0.47 


dN ch /di] 




687. 


560. 


457. 


372. 


302. 


246. 


197. 


156. 


124. 


95.3 


70.9 


52.2 


37.5 


25.6 


stat. 


error 


0.7 


0.6 


0.5 


0.5 


0.4 


0.4 


0.3 


0.3 


0.2 


0.2 


0.2 


0.1 


0.1 


0.1 


bending syst. 


error 


25. 


17. 


14. 


11. 


10. 


9.9 


9.4 


9.0 


8.2 


7.8 


7.1 


6.1 


5.2 


4.4 


full syst. 


error 


37. 


28. 


22. 


18. 


16. 


14. 


12. 


11. 


9.6 


8.6 


7.6 


6.5 


5.4 


4.5 


dN ch /dr}/(0.bN p ) 




3.89 


3.73 


3.59 


3.45 


3.34 


3.25 


3.15 


3.05 


2.96 


2.86 


2.70 


2.60 


2.48 


2.33 


bending syst. 


error 


0.14 


0.12 


0.11 


0.10 


0.12 


0.14 


0.16 


0.19 


0.21 


0.25 


0.30 


0.35 


0.41 


0.49 


full syst. 


error 


0.23 


0.22 


0.21 


0.21 


0.21 


0.22 


0.24 


0.26 


0.28 


0.32 


0.36 


0.41 


0.46 


0.55 


E T /N ch [GcV] 




0.881 


0.879 


0.881 


0.882 


0.881 


0.880 


0.875 


0.874 


0.866 


0.858 


0.851 


0.840 


0.828 


0.823 


bending syst. 


error 


0.032 


0.026 


0.021 


0.017 


0.015 


0.012 


0.011 


0.010 


0.011 


0.013 


0.017 


0.023 


0.032 


0.047 


full syst. 


error 


0.071 


0.069 


0.067 


0.066 


0.065 


0.065 


0.064 


0.064 


0.064 


0.064 


0.064 


0.065 


0.068 


0.076 



TABLE XIV: Results of the measurements by PHENIX at ^Js NN =130 GeV. Errors have the same dimension as preceding 
value. Results are plotted in Figs. and |S1 



bin [%] 




0-5 


5-10 


10-15 


15-20 


20-25 


25-30 


30-35 


35-40 


40-45 


45-50 


50-55 


55-60 


60-65 


65-70 


N p 




348. 


294. 


250. 


211. 


179. 


150. 


125. 


103. 


83.2 


66.3 


52.1 


40.1 


30.1 


21.9 


syst. 


error 


10.0 


8.9 


8.0 


7.2 


6.6 


6.0 


5.5 


5.1 


4.7 


4.3 


4.0 


3.8 


3.6 


3.4 


A ± [(m 2 ] 




138. 


123. 


110. 


99.5 


89.4 


80.6 


72.8 


65.8 


59.5 


54.3 


49.0 


44.8 


40.9 


37.4 


syst. 


error 


11. 


9.9 


8.9 


8.1 


7.3 


6.6 


6.1 


5.6 


5.2 


4.8 


4.6 


4.5 


4.6 


4.8 


dE T /dr] [GcVJ 




523. 


425. 


349. 


287. 


237. 


191. 


154. 


122. 


96.0 


73.3 


55.5 


41.0 


30.2 


21.4 


stat. 


error 


0.9 


0.9 


0.8 


0.7 


0.6 


0.0 


0.5 


0.4 


0.4 


0.3 


0.3 


0.2 


0.2 


0.2 


bending syst. 


error 


2.6 


4.2 


5.6 


7.0 


7.5 


7.6 


7.5 


7.0 


7.3 


6.2 


5.8 


5.1 


4.4 


4.1 


full syst. 


error 


27. 


22. 


19. 


16. 


14. 


12. 


11. 


9.4 


8.8 


7.3 


6.5 


5.5 


4.7 


4.2 


e Bj T [GcV fm~ 2 . 




4.7 


4.3 


3.9 


3.6 


3.3 


3.0 


2.6 


2.3 


2.0 


1.7 


1.4 


1.1 


0.9 


0.7 


full syst. 


error 


0.5 


0.5 


0.4 


0.4 


0.4 


0.3 


0.3 


0.3 


0.3 


0.2 


0.2 


0.2 


0.2 


0.2 


dE T fdr]/(0.bN p ) [GcV] 


3.01 


2.89 


2.80 


2.72 


2.65 


2.56 


2.47 


2.37 


2.31 


2.21 


2.13 


2.05 


2.01 


1.95 


bending syst. 


error 


0.01 


0.03 


0.05 


0.07 


0.09 


0.10 


0.12 


0.14 


0.19 


0.20 


0.24 


0.29 


0.34 


0.45 


full syst. 


error 


0.18 


0.18 


0.18 


0.18 


0.19 


0.20 


0.21 


0.22 


0.25 


0.26 


0.30 


0.34 


0.39 


0.49 


dN ch /di] 




602. 


488. 


403. 


329. 


270. 


219. 


176. 


139. 


109. 


84.1 


64.3 


48.4 


35.2 


25.3 


stat. 


error 


1.4 


1.2 


1.1 


0.9 


0.8 


0.7 


0.6 


0.5 


0.4 


0.3 


0.3 


0.3 


0.2 


0.2 


bending syst. 


error 


19. 


13. 


10. 


9.9 


8.6 


8.4 


8.3 


7.7 


7.5 


0.4 


5.9 


5.2 


4.3 


4.1 


full syst. 


error 


28. 


22. 


17. 


15. 


13. 


11. 


10. 


9.1 


8.4 


7.0 


6.3 


5.4 


4.5 


4.1 


dN ch /dr}/(0.5N p ) 




3.46 


3.32 


3.23 


3.12 


3.03 


2.93 


2.82 


2.70 


2.63 


2.54 


2.47 


2.41 


2.34 


2.31 


bending syst. 


error 


0.11 


0.09 


0.08 


0.09 


0.10 


0.12 


0.14 


0.16 


0.19 


0.21 


0.25 


0.30 


0.35 


0.47 


full syst. 


error 


0.19 


0.18 


0.17 


0.18 


0.18 


0.19 


0.21 


0.22 


0.25 


0.27 


0.31 


0.35 


0.41 


0.52 


Et/N cH [GeV] 




0.869 


0.870 


0.867 


0.874 


0.877 


0.873 


0.875 


0.876 


0.878 


0.871 


0.864 


0.847 


0.857 


0.844 


bending syst. 


error 


0.028 


0.023 


0.019 


0.016 


0.015 


0.014 


0.014 


0.016 


0.020 


0.025 


0.033 


0.043 


0.060 


0.083 


full syst. 


error 


0.066 


0.064 


0.063 


0.062 


0.062 


0.062 


0.062 


0.063 


0.064 


0.065 


0.068 


0.073 


0.084 


0.101 
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TABLE XV: Results of the measurements by PHENIX at ^/s NN =19.6 GeV. Errors have the same dimension as preceding 
value. Results are plotted in Figs. |S] |7| and |H] 



bin [%] 


0-5 


5-10 


10-15 


15-20 


20-25 


25-30 


30-35 


35-40 


40-45 


45-50 


N p 


336. 


288. 


243. 


204. 


172. 


144. 


120. 


98.4 


79.8 


63.8 


syst. error 


9.7 


8.8 


7.9 


7.1 


6.4 


5.9 


5.4 


5.0 


4.6 


4.3 


A , [fm 2 ! 


133 


119 


106 


95 6 


85 8 


77 2 


69 7 


62 7 


56 7 


51 3 


syst. error 


11 


g 


8 6 


7.8 


7.0 


6 4 


5.8 


5 3 


4.9 


4 6 


✓vcv, / r i, n irvvi 


230 


194 


164 


134 


109 


88 4 


72 


58 1 


45 3 


35 2 


sta/t. error 


0.8 


0.8 


0.8 


0.7 


0.5 


0.5 


0.4 


0.4 


0.3 


0.3 


bending syst error 


1.7 


2.6 


2.9 


4.0 


3.8 


3.8 


3 9 


3.9 


3.9 


3.4 


full syst . ci ror 


14 


12 


1 1 


9 3 


7.8 


6.7 


6 


5 3 


4.8 


4 


tfjj t \\j<j v mi c 


2 2 


2 


1 9 


1 8 


1 6 


1 4 


1 3 


1 2 


1 


9 


full syst. error 


0.2 


0.2 


0.2 


0.2 


0.2 


0.2 


0.2 


0.2 


0.1 


0. 1 




1.37 


1.35 


1.35 


1.31 


1.27 


1.22 


1.20 


1.18 


1.13 


1.10 


bending syst. error 


01 


02 


02 


04 


05 


05 


07 


08 


10 


11 


full syst. error 


0.09 


0.09 


0.10 


0.10 


0.10 


0.10 


0.11 


0.12 


0.14 


0.15 


dN ch /drj 


312. 


265. 


226. 


187. 


154. 


125. 


102. 


82.6 


65.0 


51.1 


stat. error 


0.9 


0.9 


0.9 


0.8 


0.7 


0.6 


0.5 


0.4 


0.4 


0.3 


bending syst. error 


5.3 


4.5 


4.4 


5.4 


4.9 


5.3 


5.2 


5.1 


5.1 


4.2 


full syst. error 


21. 


18. 


15. 


13. 


11. 


9.7 


8.5 


7.4 


6.6 


5.4 


dN ch /drif(0.5N p ) 


1.86 


1.84 


1.85 


1.83 


1.78 


1.74 


1.71 


1.68 


1.63 


1.60 


bending syst. error 


0.03 


0.03 


0.04 


0.05 


0.06 


0.08 


0.09 


0.11 


0.14 


0.14 


full syst. error 


0.14 


0.14 


0.14 


0.14 


0.14 


0.15 


0.16 


0.17 


0.19 


0.20 


E T /N ah [GeVj 


0.738 


0.733 


0.728 


0.720 


0.711 


0.705 


0.704 


0.704 


0.697 


0.690 


bending syst. error 


0.027 


0.023 


0.020 


0.017 


0.015 


0.014 


0.014 


0.016 


0.019 


0.024 


full syst. error 


0.078 


0.076 


0.075 


0.073 


0.072 


0.071 


0.071 


0.071 


0.072 


0.072 



TABLE XVI: Ratios of measured quantities at 200 GeV/130 GeV and 200 GeV/19.6 GeV. The number of N p is the average 
between two energies. The data is plotted in Fig. 



bin [%] 


0-5 


5-10 


10-15 


15-20 


20-25 


25-30 


30-35 


35-40 


40-45 


45-50 


50-55 


55-60 


60-65 


65-70 


200 GeV/130 GeV 


N p 


350. 


297. 


252. 


213. 


180. 


150. 


125. 


103. 


83.2 


66.5 


52.3 


40.2 


30.2 


22.0 


syst. error 


10. 


9.0 


8.1 


7.3 


6.6 


6.0 


5.5 


5.1 


4.7 


4.3 


4.1 


3.8 


3.6 


3.4 


dE T /dq/(0.5N p ) 


1.14 


1.13 


1.13 


1.12 


1.11 


1.12 


1.12 


1.12 


1.11 


1.11 


1.08 


1.07 


1.03 


0.98 


bending syst. error 


0.01 


0.01 


0.01 


0.01 


0.01 


0.02 


0.02 


0.03 


0.04 


0.05 


0.05 


0.07 


0.09 


0.13 


full syst. error 


0.04 


0.04 


0.04 


0.04 


0.04 


0.05 


0.05 


0.05 


0.06 


0.06 


0.07 


0.08 


0.10 


0.13 


dN ch /d v /(0.5N p ) 


1.12 


1.12 


1.11 


1.11 


1.10 


1.11 


1.12 


1.13 


1.13 


1.13 


1.10 


1.08 


1.06 


1.01 


bending syst. error 


0.01 


0.01 


0.01 


0.01 


0.01 


0.01 


0.02 


0.02 


0.03 


0.03 


0.03 


0.03 


0.04 


0.07 


full syst. error 


0.04 


0.04 


0.04 


0.04 


0.04 


0.04 


0.04 


0.04 


0.04 


0.04 


0.04 


0.05 


0.05 


0.08 


200 GeV/19.6 GeV 


N p 


344. 


294. 


249. 


210. 


177. 


148. 


122. 


101. 


81.6 


65.2 










syst. error 


10. 


9.0 


8.1 


7.3 


6.6 


6.0 


5.5 


5.1 


4.7 


4.3 










dE T /dr)/(0.5N p ) 


2.50 


2.43 


2.34 


2.32 


2.32 


2.34 


2.29 


2.25 


2.26 


2.22 










bending syst. error 


0.02 


0.03 


0.04 


0.06 


0.07 


0.09 


0.11 


0.13 


0.17 


0.19 










full syst. error 


0.09 


0.09 


0.09 


0.10 


0.10 


0.12 


0.13 


0.15 


0.19 


0.20 










dN ch /d-q/(0.5N p ) 


2.09 


2.03 


1.94 


1.89 


1.87 


1.87 


1.84 


1.81 


1.82 


1.79 










bending syst. error 


0.08 


0.06 


0.06 


0.06 


0.06 


0.07 


0.08 


0.10 


0.12 


0.13 










full syst. error 


0.15 


0.14 


0.13 


0.13 


0.13 


0.14 


0.14 


0.15 


0.17 


0.17 











TABLE XVII: Average values of dN c h/dr/ /(0.5N P ) at different ^s NN . An additional 5% error should be added to columns 
17.2 GeV through 4.8 GeV for the uncertainty related to recalculation to the Center of Mass system. The results are presented 
in Figs. 131 and [HI 





N p 


375 


350 


325 


300 


275 


250 


225 


200 


175 


150 


125 


100 


75 


50 


25 



200 GeV 


3.92 


3 


.81 


3.72 


3 


.65 


3.56 


3.51 


3.45 


3.38 


3.34 


3.27 


3.20 


3.14 


3.03 


2 


.73 


2.78 


error 


0.13 





.13 


0.12 





.12 


0.12 


0.12 


0.12 


0.11 


0.12 


0.12 


0.12 


0.13 


0.13 





.13 


0.43 


130 GeV 


3.41 


3 


.31 


3.22 


3 


.16 


3.11 


3.07 


3.04 


3.00 


2.96 


2.89 


2.83 


2.73 


2.65 


2 


.53 


2.36 


error 


0.10 





.10 


0.10 





.10 


0.09 


0.09 


0.10 


0.09 


0.10 


0.10 


0.10 


0.11 


0.11 





.12 


0.30 


19.6 GeV 




1 


.91 


1.89 


1 


.88 


1.87 


1.87 


1.85 


1.83 


1.81 


1.76 


1.72 


1.68 


1.62 








error 







.11 


0.11 





.11 


0.11 


0.12 


0.12 


0.12 


0.12 


0.13 


0.14 


0.15 


0.19 








17.2 GeV 


1.97 


1 


.93 


1.90 


1 


.88 


1.83 


1.80 


1.78 


1.75 


1.72 


1.69 


1.66 


1.66 


1.61 


1 


.54 


1.45 


error 


0.12 





.12 


0.14 





.15 


0.16 


0.17 


0.17 


0.17 


0.17 


0.17 


0.16 


0.23 


0.21 





.19 


0.13 


8.7 GeV 


1.26 


1 


.22 


1.20 


1 


.18 


1.17 


1.16 


1.16 


1.14 


1.14 


1.14 


1.14 


1.14 










error 


0.11 





.11 


0.11 





.10 


0.10 


0.10 


0.10 


0.10 


0.09 


0.09 


0.09 


0.09 










4.8 GeV 


0.92 





.89 


0.85 





.81 


0.78 


0.76 


0.75 


0.74 


0.72 


0.71 


0.70 


0.67 


0.64 





.63 




error 


0.14 





.13 


0.13 





.12 


0.12 


0.11 


0.11 


0.11 


0.11 


0.11 


0.11 


0.14 


0.18 





.21 





20 



[1] H. Hahn et al. , Nucl. Instrum. Methods A499, 245-263 
(2003). 

[2] PHENIX Collaboration, K. Adcox et al. , Nucl. Instrum. 

Methods A499, 469-479 (2003). 
[3] K. Adcox et al. , Nucl. Instrum. Methods A497, 263-293 

(2003). 

[4] L. Aphechetche et al. , Nucl. Instrum. Methods A499, 

521-536 (2003). 
[5] M. Allen et al. , Nucl. Instrum. Methods A499, 549-559 

(2003). 

[6] S. Adler et al. , Nucl. Instrum. Methods A470, 488-499 
(2001). 

[7] PHENIX Collaboration, K. Adcox et al. , Phys. Rev. 

Lett. 87, 052301 (2001). 
[8] PHENIX Collaboration, K. Adcox et al. , Phys. Rev. 

Lett. 86, 3500 (2001). 
[9] A. Milov for the PHENIX Collaboration, Nucl. 
Phys. A698, 171 (2002). 
[10] A. Bazilevsky for the PHENIX Collaboration, Nucl. 

Phys. A715, 486 (2003). 
[11] D. Silvermyr, Ph.D. Thesis. Lund University (2001). 
[12] A. Milov, Ph.D. Thesis. The Weizmann Institute of Sci- 
ence (2002). 

[13] PHENIX Collaboration, K. Adcox et al. , Phys. Rev. 

Lett. 88, 242301 (2002). 
[14] GEANT 3.2.1, CERN program library. 
[15] X.N. Wang and M. Gyulassy, Phys. Rev. D44, 3501 

(1991). 

[16] PHENIX Collaboration, S.S. Adler et al. , Phys. 

Rev. C69, 034909 (2004). 
[17] J. Bachler for the NA49 collaboration, Nucl. Phys. A661, 

45 (1999). 

[18] M. van Leeuwen for the NA49 collaboration, Nucl. 

Phys. A715, 161 (2003). 
[19] NA49 Collaboration, S.V. Afanasiev et al. , Phys. 

Rev. C66, 054902 ( 2002) 

[20] K. Hagiwara et al. , |http://pdg.lbl.gov| Phys. Rev. D66, 

010001 (2002). 

[21] B. Hahn, D.G. Ravenhall and R. Hofstadter, Phys. Rev, 
101, 1131 (1956), and 

C.W. De Jager et al. , Atomic Data and Nuclear Table 
24, 479 (1974). 

[22] PHOBOS Collaboration, B.B. Back et al. , 

nucl-ex/0210015 
[23] E-802 Collaboration, T. Abbott et al. , Phys. Rev. C63, 

064602 (2001). 

[24] PHENIX Collaboration, K. Adcox et al. , Phys. 

Rev. C66, 024901 (2002). 
[25] WA98 Collaboration, M.M. Aggarval et al. , Eur. Phys. 

J. C18, 651 (2001). 
[26] PHENIX Collaboration, S.S. Adler et al. , Phys. Rev. 

Lett. 91, 241803-1 (2003). 
[27] M. Gyulassy, pucl-th/010 6072] 

[28] A. Krasnitz, Y. Nara, R. Venugopalan, |hep-ph/0305112| 

Nucl. Phys. A727, 427-436 (2003). 
[29] J.D. Bjorken, Phys. Rev. D27, 140 (1983). 
[30] NA49 Collaboration, T. Alber et al. , Phys. Rev. Lett. 

75, 3814 (1995). 
[31] STAR Collaboration, J. Adams et al. , nucl-ex/0311017 
[32] BRAHMS Collaboration, I.G. Bearden et al. , Phys. 

Lett. B523, 227 (2001). 



[33] BRAHMS Collaboration, I.G. Bearden et al. , Phys. Rev. 

Lett. 88, 202301 (2002). 
[34] PHOBOS Collaboration, B.B. Back et al. , Phys. Rev. 

Lett. 85, 3100 (2000). 
[35] PHOBOS Collaboration, B.B. Back et al. , Phys. 

Rev. C65, 061901 (2002). 
[36] M.D. Baker et al. for the PHOBOS Collaboration,, Nucl. 

Phys. A715, 65c- 74c (2003). 
[37] T. S. Ullrich for the STAR Collaboration, 

nucl-ex/0305018 
[38] T.S. Ullrich for the STAR Collaboration, Nucl. 

Phys. A715, 399 (2003). 
[39] STAR Collaboration, J Adams et al. , nucl-ex/0407003 
[40] D. Miskowiec for the CERES Collaboration, Proc. of 

CIPPQG Palaiseau, September 2001. 
[41] NA50 Collaboration, M.C. Abreu et al. , Phys. 

Lett. B530, 43 (2002). 
[42] NA57 Collaboration, F. Antinori et al. , 

nucl-ex/0406004 
[43] PHOBOS Collaboration, B.B. Back et al. , 

nucl-ex/0301017 
[44] UA5 Collaboration, G.J. Alner et al. , Z. Phys. C33, 1-6 

(1986). 

[45] W. Thome et al. , Nucl. Phys. B129, 365-389 (1977). 
[46] X.N. Wang et al. , Phys. Rev. C68, 054902 (2003). 
[47] B. Andersson et al. , Nucl. Phys. B281, 289 (2003). 
[48] K.J. Eskola et al. , Nucl. Phys. B570, 379 (2000); Phys. 

Lett. B497, 39 (2001). 
[49] D. Kharzeev and M. Nardi, Phys. Lett. B507, 121 

(2001); D. Kharzeev and E. Levin, Phys. Lett. B523, 

79 (2001). 

[50] S. Lee and X.N. Wang, Phys. Lett. B527, 85 (2002). 

[51] Z. Lin et al. , Phys. Rev. C64, 011902 (2001). 

[52] A. Accardi, Phys. Rev. C64, 064905 (2001). 

[53] R. Ugoccioni, Phys. Lett. B491, 253 (2000). 

[54] A. Capella et al, Phys. Repts. 236, 225 (1994). 

[55] N. Armesto Perez et al. , Phys. Lett. B527, 92 (2002); 

Eur. Phys. J. C22, 149 (2001). 

[56] D.E. Kahana and S.H. Kahana, nucl-th/0208063 

[57] S. Jeon and J. Kapusta, Phys. Rev. C63, 011901 (2001). 

[58] F. Ceretto, Ph.D. Thesis. University of Heidelberg 

(1998). 

[59] H. Appelshauser for the CERES Collaboration, Nucl. 

Phys. A698, 253 (2002). 
[60] E802 Collaboration, L. Ahle et al. , Phys. Rev. C59, 

2173 (1999). 

[61] E802 Collaboration, L. Ahle et al. , Phys. Rev. C58, 

3523-3538 (1998). 
[62] C.A. Ogilvie for the E866 and E917 Collaborations, Nucl. 

Phys. A638, 57c-68c (1998). 
[63] Y. Akiba for the E802 Collaboration, Nucl. Phys. A610, 

139 (1996). 

[64] E917 Collaboration, B.B. Back et al. , Phys. Rev. Lett. 

86, 1970-1973 (2001). 
[65] E802 Collaboration, L. Ahle et al. , Phys. Rev. C57, 

R466-R470 (1998). 
[66] E802 Collaboration, L. Ahle et al. , Phys. Lett. B476, 

1 (2000). 

[67] E802 Collaboration, L. Ahle et al. , Phys. Lett. B490, 
53 (2000). 

[68] FOPI Collaboration, W. Reisdorf et al. , Nucl. 



27 



Phys. A612, 493 (1997). [70] FOPI Collaboration, B. Hong et al. , Phys. Rev. C66, 

[69] FOPI Collaboration, D. Pelte et al. , Z. Phys. A357, 034901 (2002). 

215 (1997). 



